{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "### Compare Blondin Toy Models ###" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/tqdm/autonotebook/__init__.py:14: TqdmExperimentalWarning: Using `tqdm.autonotebook.tqdm` in notebook mode. Use `tqdm.tqdm` instead to force console mode (e.g. in jupyter console)\n", " \" (e.g. in jupyter console)\", TqdmExperimentalWarning)\n" ] } ], "source": [ "import re\n", "import base64\n", "\n", "import pandas as pd\n", "import numpy as np\n", "\n", "import yaml\n", "\n", "from tardis import run_tardis\n", "from astropy import units as u\n", "\n", "from IPython.display import HTML\n", "\n", "from tardis.io.parsers import read_blondin_toymodel, convert_blondin_toymodel\n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "\n", "def create_df_download_link( df, title=\"Download CSV file\", filename=\"data.csv\", type='zip'):\n", " if type != 'zip':\n", " csv = df.to_csv()\n", " b64 = base64.b64encode(csv.encode())\n", " payload = b64.decode()\n", " else:\n", " payload = df.decode()\n", " html = '{title}'\n", " html = html.format(payload=payload,title=title,filename=filename)\n", " return HTML(html)\n", "\n", "def create_download_link(download_data, title=\"Download CSV file\", filename=\"data.csv\", type='zip'):\n", " if type != 'zip':\n", " csv = df.to_csv()\n", " b64 = base64.b64encode(csv.encode())\n", " payload = b64.decode()\n", " else:\n", " payload = df.decode()\n", " html = '{title}'\n", " html = html.format(payload=payload,title=title,filename=filename)\n", " return HTML(html)\n", "\n", "pattern_remove_bracket = re.compile('\\[.+\\]')\n", "t0_pattern = re.compile('tend = (.+)\\n')" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "#### Comparison of snia_toy01.dat ####" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "convert_blondin_toymodel('snia_toy01.dat', 'snia_toy01.test.csvy', 9000, 30000, conversion_t_electron_rad=1/0.9)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "!code snia_toy01.test.csvy" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "blondin_dict['v_inner_boundary'] = '9000 km/s'\n", "blondin_dict['v_outer_boundary'] = '20000 km/s'" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "csvy_file = '---\\n{0}\\n---\\n{1}'.format(yaml.dump(blondin_dict, default_flow_style=False), blondin_csv.to_csv(index=False))" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "with open('blondin_compare_01.csvy', 'w') as fh:\n", " fh.write(csvy_file)" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "**Download the model and setup**\n", "\n", "[blondin_compare_01.csvy](blondin_compare_01.csvy)\n", "\n", "[blondin_model_compare_01.yml](blondin_model_compare_01.yml)\n", "\n", "**Assumptions and explanations**\n", "\n", "We assume $T_\\textrm{electron} = 0.9 T_\\textrm{radiative}$\n", "\n", "\n", "[Plasma assumptions](../../../physics/setup/plasma/nebular_plasma.rst)\n", "\n", "\n", "[Line interaction assumptions](../../../physics/montecarlo/lineinteraction.rst)\n", "\n" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/importlib/_bootstrap.py:219: QAWarning: pyne.data is not yet QA compliant.\n", " return f(*args, **kwds)\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/importlib/_bootstrap.py:219: QAWarning: pyne.material is not yet QA compliant.\n", " return f(*args, **kwds)\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.\n", " AstropyDeprecationWarning)\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mtardis.plasma.standard_plasmas\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Reading Atomic Data from kurucz_cd23_chianti_H_He.h5 (\u001b[1mstandard_plasmas.py\u001b[0m:74)\n", "[\u001b[1mtardis.io.atom_data.util\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Atom Data kurucz_cd23_chianti_H_He.h5 not found in local path. Exists in TARDIS Data repo /Users/wkerzend/projects/tardis/tardis-data/kurucz_cd23_chianti_H_He.h5 (\u001b[1mutil.py\u001b[0m:29)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/IPython/core/interactiveshell.py:3267: PerformanceWarning: indexing past lexsort depth may impact performance.\n", " exec(code_obj, self.user_global_ns, self.user_ns)\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mtardis.io.atom_data.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Read Atom Data with UUID=6f7b09e887a311e7a06b246e96350010 and MD5=864f1753714343c41f99cb065710cace. (\u001b[1mbase.py\u001b[0m:184)\n", "[\u001b[1mtardis.io.atom_data.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Non provided atomic data: synpp_refs, photoionization_data (\u001b[1mbase.py\u001b[0m:187)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.\n", " AstropyDeprecationWarning)\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/python/tardis/tardis/plasma/properties/atomic.py:195: FutureWarning: .labels was deprecated in version 0.24.0. Use .codes instead.\n", " zeta_data['atomic_number'] = zeta_data.index.labels[0] + 1\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/python/tardis/tardis/plasma/properties/atomic.py:196: FutureWarning: .labels was deprecated in version 0.24.0. Use .codes instead.\n", " zeta_data['ion_number'] = zeta_data.index.labels[1] + 1\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mtardis.plasma.properties.atomic\u001b[0m][\u001b[1;33mWARNING\u001b[0m] Zeta_data missing - replaced with 1s. Missing ions: [(14, 15), (16, 17), (20, 21), (26, 27), (27, 28), (28, 29)] (\u001b[1matomic.py\u001b[0m:215)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/python/tardis/tardis/plasma/properties/ion_population.py:63: FutureWarning: \n", "Passing list-likes to .loc or [] with any missing label will raise\n", "KeyError in the future, you can use .reindex() as an alternative.\n", "\n", "See the documentation here:\n", "https://pandas.pydata.org/pandas-docs/stable/indexing.html#deprecate-loc-reindex-listlike\n", " partition_function.index].dropna())\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/equivalencies.py:90: RuntimeWarning: divide by zero encountered in double_scalars\n", " (si.m, si.Hz, lambda x: _si.c.value / x),\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/astropy/units/quantity.py:1067: AstropyDeprecationWarning: The truth value of a Quantity is ambiguous. In the future this will raise a ValueError.\n", " AstropyDeprecationWarning)\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 1/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.30150e+43 erg / s Luminosity absorbed = 7.43372e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 5000.0 7831.059753 0.462810 0.987065\n", "\t5 5000.0 7106.913190 0.379149 1.304910\n", "\t10 5000.0 6779.027752 0.336294 1.473708\n", "\t15 5000.0 6591.160845 0.304883 1.556553\n", "\t20 5000.0 6469.413444 0.279753 1.549488\n", "\t25 5000.0 6369.832186 0.258765 1.548078\n", "\t30 5000.0 6285.178795 0.240774 1.531938\n", "\t35 5000.0 6216.656067 0.225078 1.509501\n", "\t40 5000.0 6155.241379 0.211205 1.474743\n", "\t45 5000.0 6105.918970 0.198821 1.414270\n", "\t50 5000.0 6076.371342 0.187679 1.348556\n", "\t55 5000.0 6030.241880 0.177588 1.311088\n", "\t60 5000.0 6005.973722 0.168399 1.256169\n", "\t65 5000.0 6002.125995 0.159992 1.169288\n", "\t70 5000.0 5976.207283 0.152269 1.122986\n", "\t75 5000.0 5951.666260 0.145148 1.077897\n", "\t80 5000.0 5918.254016 0.138563 1.053437\n", "\t85 5000.0 5918.208740 0.132453 0.989942\n", "\t90 5000.0 5867.100159 0.126772 0.990185\n", "\t95 5000.0 5859.007030 0.121475 0.929519\n", "\t100 5000.0 5843.254166 0.116527 0.894981\n", "\t105 5000.0 5841.131195 0.111894 0.849290\n", "\t110 5000.0 5821.292215 0.107550 0.811782\n", "\t115 5000.0 5820.272035 0.103468 0.769661\n", "\t120 5000.0 5811.350827 0.099626 0.736181\n", "\t125 5000.0 5793.908583 0.096006 0.716442\n", "\t130 5000.0 5795.478604 0.092589 0.681608\n", "\t135 5000.0 5773.849040 0.089359 0.667536\n", "\t140 5000.0 5775.403713 0.086303 0.634047\n", "\t145 5000.0 5736.539129 0.083408 0.630179\n", "\t150 5000.0 5732.357290 0.080662 0.601033\n", "\t155 5000.0 5743.364284 0.078054 0.570580\n", "\t160 5000.0 5733.145769 0.075575 0.551961\n", "\t165 5000.0 5718.331258 0.073217 0.529948\n", "\t170 5000.0 5699.632837 0.070971 0.519324\n", "\t175 5000.0 5695.349689 0.068830 0.499699\n", "\t180 5000.0 5688.534677 0.066787 0.478921\n", "\t185 5000.0 5683.055043 0.064837 0.463320\n", "\t190 5000.0 5665.611101 0.062973 0.450357\n", "\t195 5000.0 5654.224456 0.061191 0.437244\n", "\t200 5000.0 5651.225192 0.059485 0.421169\n", "\t205 5000.0 5626.660446 0.057852 0.414763\n", "\t210 5000.0 5595.199895 0.056286 0.409042\n", "\t215 5000.0 5577.506895 0.054785 0.400196\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/python/tardis/tardis/plasma/properties/ion_population.py:63: FutureWarning: \n", "Passing list-likes to .loc or [] with any missing label will raise\n", "KeyError in the future, you can use .reindex() as an alternative.\n", "\n", "See the documentation here:\n", "https://pandas.pydata.org/pandas-docs/stable/indexing.html#deprecate-loc-reindex-listlike\n", " partition_function.index].dropna())\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 2/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.38526e+43 erg / s Luminosity absorbed = 6.66321e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 7831.059753 9943.504253 0.987065 0.941908\n", "\t5 7106.913190 9168.067423 1.304910 1.150356\n", "\t10 6779.027752 8678.999816 1.473708 1.286543\n", "\t15 6591.160845 8377.317578 1.556553 1.346125\n", "\t20 6469.413444 8151.787425 1.549488 1.358795\n", "\t25 6369.832186 7985.142325 1.548078 1.353925\n", "\t30 6285.178795 7878.847609 1.531938 1.319217\n", "\t35 6216.656067 7741.589421 1.509501 1.307884\n", "\t40 6155.241379 7655.316512 1.474743 1.270084\n", "\t45 6105.918970 7556.986602 1.414270 1.234470\n", "\t50 6076.371342 7522.785118 1.348556 1.170937\n", "\t55 6030.241880 7447.391749 1.311088 1.133516\n", "\t60 6005.973722 7408.925630 1.256169 1.089917\n", "\t65 6002.125995 7364.361995 1.169288 1.029933\n", "\t70 5976.207283 7336.100563 1.122986 0.990663\n", "\t75 5951.666260 7306.592489 1.077897 0.946869\n", "\t80 5918.254016 7238.893988 1.053437 0.931127\n", "\t85 5918.208740 7226.363390 0.989942 0.877828\n", "\t90 5867.100159 7197.898366 0.990185 0.857879\n", "\t95 5859.007030 7165.997446 0.929519 0.812933\n", "\t100 5843.254166 7159.070470 0.894981 0.775213\n", "\t105 5841.131195 7143.741089 0.849290 0.738824\n", "\t110 5821.292215 7119.704519 0.811782 0.707784\n", "\t115 5820.272035 7112.481717 0.769661 0.671211\n", "\t120 5811.350827 7060.505651 0.736181 0.652498\n", "\t125 5793.908583 7066.008522 0.716442 0.624559\n", "\t130 5795.478604 7029.029121 0.681608 0.605209\n", "\t135 5773.849040 7021.659449 0.667536 0.582533\n", "\t140 5775.403713 7010.212505 0.634047 0.555628\n", "\t145 5736.539129 6990.562369 0.630179 0.542516\n", "\t150 5732.357290 6984.884673 0.601033 0.518933\n", "\t155 5743.364284 6980.471223 0.570580 0.494955\n", "\t160 5733.145769 6938.386837 0.551961 0.484650\n", "\t165 5718.331258 6928.177798 0.529948 0.465554\n", "\t170 5699.632837 6904.323015 0.519324 0.451990\n", "\t175 5695.349689 6905.684428 0.499699 0.434141\n", "\t180 5688.534677 6893.408541 0.478921 0.418797\n", "\t185 5683.055043 6881.358878 0.463320 0.404739\n", "\t190 5665.611101 6848.434616 0.450357 0.394703\n", "\t195 5654.224456 6831.751052 0.437244 0.382167\n", "\t200 5651.225192 6821.557616 0.421169 0.367724\n", "\t205 5626.660446 6798.252052 0.414763 0.360225\n", "\t210 5595.199895 6760.804314 0.409042 0.354866\n", "\t215 5577.506895 6740.013534 0.400196 0.344957\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 3/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.32445e+43 erg / s Luminosity absorbed = 7.18316e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 9943.504253 11184.865494 0.941908 0.874243\n", "\t5 9168.067423 10590.890516 1.150356 0.967498\n", "\t10 8678.999816 10184.538479 1.286543 1.029334\n", "\t15 8377.317578 9900.090529 1.346125 1.054300\n", "\t20 8151.787425 9664.894816 1.358795 1.054968\n", "\t25 7985.142325 9510.226327 1.353925 1.037032\n", "\t30 7878.847609 9394.106182 1.319217 1.004565\n", "\t35 7741.589421 9271.496145 1.307884 0.982576\n", "\t40 7655.316512 9161.047735 1.270084 0.949923\n", "\t45 7556.986602 9072.338911 1.234470 0.915208\n", "\t50 7522.785118 9010.672618 1.170937 0.872802\n", "\t55 7447.391749 8901.383373 1.133516 0.845834\n", "\t60 7408.925630 8834.036103 1.089917 0.815912\n", "\t65 7364.361995 8769.307119 1.029933 0.772520\n", "\t70 7336.100563 8711.903127 0.990663 0.746307\n", "\t75 7306.592489 8640.885761 0.946869 0.718271\n", "\t80 7238.893988 8542.892115 0.931127 0.704924\n", "\t85 7226.363390 8516.723270 0.877828 0.666870\n", "\t90 7197.898366 8502.592440 0.857879 0.642794\n", "\t95 7165.997446 8447.746462 0.812933 0.611320\n", "\t100 7159.070470 8412.121792 0.775213 0.587311\n", "\t105 7143.741089 8381.138885 0.738824 0.562031\n", "\t110 7119.704519 8359.007818 0.707784 0.537153\n", "\t115 7112.481717 8322.937242 0.671211 0.513156\n", "\t120 7060.505651 8249.430763 0.652498 0.500380\n", "\t125 7066.008522 8275.891240 0.624559 0.473512\n", "\t130 7029.029121 8204.492926 0.605209 0.461398\n", "\t135 7021.659449 8180.726449 0.582533 0.445042\n", "\t140 7010.212505 8169.127793 0.555628 0.424694\n", "\t145 6990.562369 8158.252968 0.542516 0.411711\n", "\t150 6984.884673 8112.858348 0.518933 0.397343\n", "\t155 6980.471223 8098.797522 0.494955 0.379947\n", "\t160 6938.386837 8055.533873 0.484650 0.370261\n", "\t165 6928.177798 8055.402146 0.465554 0.354906\n", "\t170 6904.323015 8019.347064 0.451990 0.345449\n", "\t175 6905.684428 8006.258863 0.434141 0.331780\n", "\t180 6893.408541 8002.951619 0.418797 0.319606\n", "\t185 6881.358878 7970.659310 0.404739 0.309280\n", "\t190 6848.434616 7932.276565 0.394703 0.301617\n", "\t195 6831.751052 7915.200790 0.382167 0.291916\n", "\t200 6821.557616 7898.840956 0.367724 0.280425\n", "\t205 6798.252052 7883.486836 0.360225 0.273069\n", "\t210 6760.804314 7840.128677 0.354866 0.268436\n", "\t215 6740.013534 7815.699388 0.344957 0.260891\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 4/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.30498e+43 erg / s Luminosity absorbed = 7.38203e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 11184.865494 11911.711939 0.874243 0.815342\n", "\t5 10590.890516 11423.170586 0.967498 0.852032\n", "\t10 10184.538479 11134.470932 1.029334 0.863305\n", "\t15 9900.090529 10887.349224 1.054300 0.863750\n", "\t20 9664.894816 10740.461976 1.054968 0.832290\n", "\t25 9510.226327 10562.641071 1.037032 0.813050\n", "\t30 9394.106182 10481.388048 1.004565 0.774264\n", "\t35 9271.496145 10372.742105 0.982576 0.747951\n", "\t40 9161.047735 10274.337236 0.949923 0.716987\n", "\t45 9072.338911 10202.946659 0.915208 0.684430\n", "\t50 9010.672618 10120.383319 0.872802 0.652331\n", "\t55 8901.383373 10039.433120 0.845834 0.624391\n", "\t60 8834.036103 9985.105525 0.815912 0.597331\n", "\t65 8769.307119 9909.687407 0.772520 0.567351\n", "\t70 8711.903127 9853.157362 0.746307 0.545114\n", "\t75 8640.885761 9803.411756 0.718271 0.521662\n", "\t80 8542.892115 9703.094119 0.704924 0.508055\n", "\t85 8516.723270 9664.369661 0.666870 0.483056\n", "\t90 8502.592440 9644.808583 0.642794 0.462651\n", "\t95 8447.746462 9595.659098 0.611320 0.439338\n", "\t100 8412.121792 9572.851164 0.587311 0.421222\n", "\t105 8381.138885 9519.270385 0.562031 0.403866\n", "\t110 8359.007818 9520.334232 0.537153 0.383152\n", "\t115 8322.937242 9482.738472 0.513156 0.366423\n", "\t120 8249.430763 9426.321672 0.500380 0.356164\n", "\t125 8275.891240 9391.728992 0.473512 0.340621\n", "\t130 8204.492926 9367.333274 0.461398 0.327122\n", "\t135 8180.726449 9303.898069 0.445042 0.318214\n", "\t140 8169.127793 9283.891027 0.424694 0.304527\n", "\t145 8158.252968 9270.457341 0.411711 0.293142\n", "\t150 8112.858348 9230.938730 0.397343 0.282674\n", "\t155 8098.797522 9190.155128 0.379947 0.272658\n", "\t160 8055.533873 9192.939613 0.370261 0.262679\n", "\t165 8055.402146 9166.640157 0.354906 0.253095\n", "\t170 8019.347064 9116.658009 0.345449 0.246487\n", "\t175 8006.258863 9096.261076 0.331780 0.237289\n", "\t180 8002.951619 9095.824747 0.319606 0.228009\n", "\t185 7970.659310 9083.284785 0.309280 0.219454\n", "\t190 7932.276565 9051.094660 0.301617 0.212873\n", "\t195 7915.200790 9020.044932 0.291916 0.207175\n", "\t200 7898.840956 8998.843933 0.280425 0.199036\n", "\t205 7883.486836 8982.946944 0.273069 0.193372\n", "\t210 7840.128677 8952.622269 0.268436 0.188983\n", "\t215 7815.699388 8935.684680 0.260891 0.182965\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 5/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.28254e+43 erg / s Luminosity absorbed = 7.59247e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 11911.711939 12274.410751 0.815342 0.789790\n", "\t5 11423.170586 11871.455425 0.852032 0.792676\n", "\t10 11134.470932 11649.695187 0.863305 0.773668\n", "\t15 10887.349224 11471.210556 0.863750 0.748405\n", "\t20 10740.461976 11360.332943 0.832290 0.712461\n", "\t25 10562.641071 11218.786535 0.813050 0.682594\n", "\t30 10481.388048 11129.781624 0.774264 0.647152\n", "\t35 10372.742105 11054.982577 0.747951 0.615252\n", "\t40 10274.337236 10978.674047 0.716987 0.581735\n", "\t45 10202.946659 10896.716952 0.684430 0.554543\n", "\t50 10120.383319 10820.747796 0.652331 0.526362\n", "\t55 10039.433120 10788.603828 0.624391 0.494951\n", "\t60 9985.105525 10721.956293 0.597331 0.471665\n", "\t65 9909.687407 10655.755485 0.567351 0.446374\n", "\t70 9853.157362 10634.209059 0.545114 0.423620\n", "\t75 9803.411756 10590.884258 0.521662 0.402309\n", "\t80 9703.094119 10508.509946 0.508055 0.388143\n", "\t85 9664.369661 10508.882030 0.483056 0.365331\n", "\t90 9644.808583 10483.115608 0.462651 0.347868\n", "\t95 9595.659098 10430.871132 0.439338 0.329624\n", "\t100 9572.851164 10415.185286 0.421222 0.314476\n", "\t105 9519.270385 10368.229649 0.403866 0.300532\n", "\t110 9520.334232 10371.842110 0.383152 0.284290\n", "\t115 9482.738472 10344.461437 0.366423 0.270853\n", "\t120 9426.321672 10330.632364 0.356164 0.260494\n", "\t125 9391.728992 10280.580995 0.340621 0.250351\n", "\t130 9367.333274 10249.596303 0.327122 0.239675\n", "\t135 9303.898069 10236.391898 0.318214 0.229879\n", "\t140 9283.891027 10200.349641 0.304527 0.220010\n", "\t145 9270.457341 10180.401007 0.293142 0.212525\n", "\t150 9230.938730 10140.384246 0.282674 0.204939\n", "\t155 9190.155128 10095.159566 0.272658 0.197206\n", "\t160 9192.939613 10116.202875 0.262679 0.188372\n", "\t165 9166.640157 10064.971970 0.253095 0.182719\n", "\t170 9116.658009 10059.156300 0.246487 0.176322\n", "\t175 9096.261076 10026.674915 0.237289 0.170347\n", "\t180 9095.824747 9986.590882 0.228009 0.165054\n", "\t185 9083.284785 9991.204987 0.219454 0.157743\n", "\t190 9051.094660 9954.915111 0.212873 0.153163\n", "\t195 9020.044932 9946.636259 0.207175 0.148020\n", "\t200 8998.843933 9904.779162 0.199036 0.142722\n", "\t205 8982.946944 9878.309679 0.193372 0.138820\n", "\t210 8952.622269 9865.404965 0.188983 0.134775\n", "\t215 8935.684680 9838.418760 0.182965 0.130692\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 6/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.28283e+43 erg / s Luminosity absorbed = 7.59368e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 12274.410751 12621.985366 0.789790 0.740187\n", "\t5 11871.455425 12288.282401 0.792676 0.714197\n", "\t10 11649.695187 12107.496695 0.773668 0.683174\n", "\t15 11471.210556 11966.421542 0.748405 0.651469\n", "\t20 11360.332943 11869.120463 0.712461 0.610504\n", "\t25 11218.786535 11782.584290 0.682594 0.576377\n", "\t30 11129.781624 11661.517833 0.647152 0.549776\n", "\t35 11054.982577 11596.680082 0.615252 0.515964\n", "\t40 10978.674047 11486.100778 0.581735 0.492717\n", "\t45 10896.716952 11411.634100 0.554543 0.466075\n", "\t50 10820.747796 11320.929647 0.526362 0.443230\n", "\t55 10788.603828 11304.819988 0.494951 0.411733\n", "\t60 10721.956293 11223.787471 0.471665 0.391653\n", "\t65 10655.755485 11187.391580 0.446374 0.366972\n", "\t70 10634.209059 11145.379508 0.423620 0.349026\n", "\t75 10590.884258 11132.155630 0.402309 0.327809\n", "\t80 10508.509946 11080.451217 0.388143 0.312311\n", "\t85 10508.882030 11053.402190 0.365331 0.294910\n", "\t90 10483.115608 11022.793935 0.347868 0.280215\n", "\t95 10430.871132 11018.638959 0.329624 0.263074\n", "\t100 10415.185286 10969.698101 0.314476 0.251057\n", "\t105 10368.229649 10996.575238 0.300532 0.235733\n", "\t110 10371.842110 10946.225411 0.284290 0.225373\n", "\t115 10344.461437 10915.429238 0.270853 0.215369\n", "\t120 10330.632364 10911.516887 0.260494 0.204801\n", "\t125 10280.580995 10876.787083 0.250351 0.195768\n", "\t130 10249.596303 10846.880881 0.239675 0.187419\n", "\t135 10236.391898 10827.878807 0.229879 0.179927\n", "\t140 10200.349641 10843.147042 0.220010 0.169929\n", "\t145 10180.401007 10836.452058 0.212525 0.163079\n", "\t150 10140.384246 10779.508126 0.204939 0.157783\n", "\t155 10095.159566 10792.246053 0.197206 0.150208\n", "\t160 10116.202875 10777.101313 0.188372 0.144407\n", "\t165 10064.971970 10764.279483 0.182719 0.138515\n", "\t170 10059.156300 10755.323524 0.176322 0.133355\n", "\t175 10026.674915 10716.093969 0.170347 0.129014\n", "\t180 9986.590882 10681.240209 0.165054 0.124735\n", "\t185 9991.204987 10710.349167 0.157743 0.118591\n", "\t190 9954.915111 10647.113250 0.153163 0.115792\n", "\t195 9946.636259 10621.085728 0.148020 0.112042\n", "\t200 9904.779162 10611.632892 0.142722 0.107582\n", "\t205 9878.309679 10575.876336 0.138820 0.104696\n", "\t210 9865.404965 10534.592000 0.134775 0.102013\n", "\t215 9838.418760 10485.222357 0.130692 0.099331\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 7/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.28515e+43 erg / s Luminosity absorbed = 7.57676e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 12621.985366 12759.629863 0.740187 0.721720\n", "\t5 12288.282401 12486.425455 0.714197 0.679950\n", "\t10 12107.496695 12369.251952 0.683174 0.636363\n", "\t15 11966.421542 12252.854693 0.651469 0.595850\n", "\t20 11869.120463 12167.826060 0.610504 0.556063\n", "\t25 11782.584290 12071.481834 0.576377 0.522488\n", "\t30 11661.517833 11955.870966 0.549776 0.496104\n", "\t35 11596.680082 11917.514987 0.515964 0.459262\n", "\t40 11486.100778 11870.688289 0.492717 0.427744\n", "\t45 11411.634100 11790.779683 0.466075 0.405292\n", "\t50 11320.929647 11703.842476 0.443230 0.384510\n", "\t55 11304.819988 11687.592974 0.411733 0.356790\n", "\t60 11223.787471 11618.853302 0.391653 0.337614\n", "\t65 11187.391580 11585.501944 0.366972 0.314979\n", "\t70 11145.379508 11526.322819 0.349026 0.299224\n", "\t75 11132.155630 11523.858835 0.327809 0.281351\n", "\t80 11080.451217 11448.221615 0.312311 0.269039\n", "\t85 11053.402190 11462.009705 0.294910 0.249884\n", "\t90 11022.793935 11416.159361 0.280215 0.238194\n", "\t95 11018.638959 11406.804183 0.263074 0.222978\n", "\t100 10969.698101 11376.225310 0.251057 0.211479\n", "\t105 10996.575238 11433.350181 0.235733 0.196343\n", "\t110 10946.225411 11372.910863 0.225373 0.189059\n", "\t115 10915.429238 11332.677175 0.215369 0.180570\n", "\t120 10911.516887 11285.604610 0.204801 0.172781\n", "\t125 10876.787083 11299.429914 0.195768 0.162189\n", "\t130 10846.880881 11263.153673 0.187419 0.155143\n", "\t135 10827.878807 11204.982403 0.179927 0.150862\n", "\t140 10843.147042 11258.798865 0.169929 0.141258\n", "\t145 10836.452058 11231.762696 0.163079 0.135271\n", "\t150 10779.508126 11207.707717 0.157783 0.129487\n", "\t155 10792.246053 11191.745451 0.150208 0.124238\n", "\t160 10777.101313 11181.509283 0.144407 0.119067\n", "\t165 10764.279483 11162.391880 0.138515 0.114391\n", "\t170 10755.323524 11158.541462 0.133355 0.109545\n", "\t175 10716.093969 11151.075575 0.129014 0.105748\n", "\t180 10681.240209 11118.631655 0.124735 0.102085\n", "\t185 10710.349167 11103.407863 0.118591 0.097843\n", "\t190 10647.113250 11072.785328 0.115792 0.094608\n", "\t195 10621.085728 11057.565986 0.112042 0.091442\n", "\t200 10611.632892 11035.369976 0.107582 0.088212\n", "\t205 10575.876336 11002.115565 0.104696 0.085741\n", "\t210 10534.592000 10955.610971 0.102013 0.083497\n", "\t215 10485.222357 10917.579162 0.099331 0.081022\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 8/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.27712e+43 erg / s Luminosity absorbed = 7.65528e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 12759.629863 12860.359189 0.721720 0.708992\n", "\t5 12486.425455 12638.631249 0.679950 0.651342\n", "\t10 12369.251952 12561.844673 0.636363 0.598991\n", "\t15 12252.854693 12457.507444 0.595850 0.554747\n", "\t20 12167.826060 12405.099383 0.556063 0.514695\n", "\t25 12071.481834 12309.733919 0.522488 0.482479\n", "\t30 11955.870966 12245.489227 0.496104 0.450183\n", "\t35 11917.514987 12183.541188 0.459262 0.418285\n", "\t40 11870.688289 12158.321534 0.427744 0.386364\n", "\t45 11790.779683 12097.953291 0.405292 0.362136\n", "\t50 11703.842476 12043.102357 0.384510 0.340574\n", "\t55 11687.592974 11992.235264 0.356790 0.317048\n", "\t60 11618.853302 11931.000829 0.337614 0.297179\n", "\t65 11585.501944 11868.817978 0.314979 0.281821\n", "\t70 11526.322819 11845.559084 0.299224 0.264268\n", "\t75 11523.858835 11839.159985 0.281351 0.247294\n", "\t80 11448.221615 11766.869353 0.269039 0.236873\n", "\t85 11462.009705 11781.216970 0.249884 0.219222\n", "\t90 11416.159361 11763.657217 0.238194 0.206710\n", "\t95 11406.804183 11736.050598 0.222978 0.194684\n", "\t100 11376.225310 11709.075282 0.211479 0.184521\n", "\t105 11433.350181 11764.053125 0.196343 0.171069\n", "\t110 11372.910863 11693.485667 0.189059 0.164673\n", "\t115 11332.677175 11662.729622 0.180570 0.157151\n", "\t120 11285.604610 11626.621706 0.172781 0.149390\n", "\t125 11299.429914 11591.005256 0.162189 0.141987\n", "\t130 11263.153673 11621.706183 0.155143 0.133331\n", "\t135 11204.982403 11562.094749 0.150862 0.129172\n", "\t140 11258.798865 11593.366562 0.141258 0.121613\n", "\t145 11231.762696 11543.895750 0.135271 0.117417\n", "\t150 11207.707717 11510.519401 0.129487 0.112442\n", "\t155 11191.745451 11521.475871 0.124238 0.107162\n", "\t160 11181.509283 11508.760099 0.119067 0.102595\n", "\t165 11162.391880 11476.537118 0.114391 0.098371\n", "\t170 11158.541462 11464.294275 0.109545 0.094355\n", "\t175 11151.075575 11445.541800 0.105748 0.091276\n", "\t180 11118.631655 11422.237088 0.102085 0.087966\n", "\t185 11103.407863 11397.864821 0.097843 0.084664\n", "\t190 11072.785328 11382.914370 0.094608 0.081583\n", "\t195 11057.565986 11349.784124 0.091442 0.078976\n", "\t200 11035.369976 11322.806026 0.088212 0.076562\n", "\t205 11002.115565 11263.555912 0.085741 0.074833\n", "\t210 10955.610971 11221.660698 0.083497 0.072828\n", "\t215 10917.579162 11200.040416 0.081022 0.070359\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 9/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.26945e+43 erg / s Luminosity absorbed = 7.72465e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 12860.359189 12957.621989 0.708992 0.691150\n", "\t5 12638.631249 12803.884316 0.651342 0.624900\n", "\t10 12561.844673 12730.751370 0.598991 0.569845\n", "\t15 12457.507444 12604.055249 0.554747 0.528149\n", "\t20 12405.099383 12583.906104 0.514695 0.485502\n", "\t25 12309.733919 12491.390514 0.482479 0.453408\n", "\t30 12245.489227 12424.403292 0.450183 0.423553\n", "\t35 12183.541188 12380.236496 0.418285 0.388916\n", "\t40 12158.321534 12363.358789 0.386364 0.359253\n", "\t45 12097.953291 12328.785784 0.362136 0.335280\n", "\t50 12043.102357 12251.735081 0.340574 0.314188\n", "\t55 11992.235264 12261.737757 0.317048 0.288563\n", "\t60 11931.000829 12152.175154 0.297179 0.274420\n", "\t65 11868.817978 12117.845990 0.281821 0.257238\n", "\t70 11845.559084 12093.081412 0.264268 0.241895\n", "\t75 11839.159985 12093.052004 0.247294 0.224445\n", "\t80 11766.869353 12037.640748 0.236873 0.214143\n", "\t85 11781.216970 12049.271412 0.219222 0.198544\n", "\t90 11763.657217 12015.545530 0.206710 0.188042\n", "\t95 11736.050598 11993.213751 0.194684 0.176258\n", "\t100 11709.075282 11930.668715 0.184521 0.168432\n", "\t105 11764.053125 12037.937411 0.171069 0.153223\n", "\t110 11693.485667 11988.734420 0.164673 0.146614\n", "\t115 11662.729622 11937.603116 0.157151 0.140315\n", "\t120 11626.621706 11887.030622 0.149390 0.133416\n", "\t125 11591.005256 11836.102537 0.141987 0.127737\n", "\t130 11621.706183 11833.294958 0.133331 0.120973\n", "\t135 11562.094749 11777.150834 0.129172 0.116478\n", "\t140 11593.366562 11825.884043 0.121613 0.109144\n", "\t145 11543.895750 11766.196357 0.117417 0.105560\n", "\t150 11510.519401 11740.623271 0.112442 0.101324\n", "\t155 11521.475871 11732.541118 0.107162 0.096637\n", "\t160 11508.760099 11708.165413 0.102595 0.092791\n", "\t165 11476.537118 11672.422129 0.098371 0.089207\n", "\t170 11464.294275 11660.830520 0.094355 0.085588\n", "\t175 11445.541800 11633.201479 0.091276 0.082922\n", "\t180 11422.237088 11615.548095 0.087966 0.079616\n", "\t185 11397.864821 11566.309428 0.084664 0.077288\n", "\t190 11382.914370 11575.749336 0.081583 0.073956\n", "\t195 11349.784124 11543.883213 0.078976 0.071674\n", "\t200 11322.806026 11495.598254 0.076562 0.069762\n", "\t205 11263.555912 11490.479524 0.074833 0.067104\n", "\t210 11221.660698 11428.432630 0.072828 0.065627\n", "\t215 11200.040416 11408.594202 0.070359 0.063437\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 10/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.26942e+43 erg / s Luminosity absorbed = 7.72233e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 12957.621989 13015.569839 0.691150 0.682918\n", "\t5 12803.884316 12904.875006 0.624900 0.607993\n", "\t10 12730.751370 12842.565429 0.569845 0.550057\n", "\t15 12604.055249 12753.343916 0.528149 0.505693\n", "\t20 12583.906104 12730.133218 0.485502 0.462479\n", "\t25 12491.390514 12664.083226 0.453408 0.427473\n", "\t30 12424.403292 12588.842740 0.423553 0.399502\n", "\t35 12380.236496 12576.400081 0.388916 0.365150\n", "\t40 12363.358789 12571.096532 0.359253 0.335182\n", "\t45 12328.785784 12567.267055 0.335280 0.308863\n", "\t50 12251.735081 12490.762479 0.314188 0.289858\n", "\t55 12261.737757 12473.100822 0.288563 0.268680\n", "\t60 12152.175154 12403.844620 0.274420 0.252799\n", "\t65 12117.845990 12379.017248 0.257238 0.235858\n", "\t70 12093.081412 12384.414417 0.241895 0.219277\n", "\t75 12093.052004 12359.661933 0.224445 0.205534\n", "\t80 12037.640748 12301.567714 0.214143 0.194798\n", "\t85 12049.271412 12282.445145 0.198544 0.182032\n", "\t90 12015.545530 12245.263730 0.188042 0.172696\n", "\t95 11993.213751 12240.211097 0.176258 0.161126\n", "\t100 11930.668715 12166.777597 0.168432 0.154560\n", "\t105 12037.937411 12230.213020 0.153223 0.141698\n", "\t110 11988.734420 12176.872696 0.146614 0.135522\n", "\t115 11937.603116 12173.728933 0.140315 0.128104\n", "\t120 11887.030622 12135.497772 0.133416 0.121988\n", "\t125 11836.102537 12064.394959 0.127737 0.117167\n", "\t130 11833.294958 12010.587301 0.120973 0.112690\n", "\t135 11777.150834 11986.931983 0.116478 0.107698\n", "\t140 11825.884043 12025.244790 0.109144 0.100791\n", "\t145 11766.196357 11995.910253 0.105560 0.096670\n", "\t150 11740.623271 11947.593265 0.101324 0.093294\n", "\t155 11732.541118 11957.542102 0.096637 0.088421\n", "\t160 11708.165413 11910.670130 0.092791 0.085384\n", "\t165 11672.422129 11888.528792 0.089207 0.081976\n", "\t170 11660.830520 11879.429183 0.085588 0.078285\n", "\t175 11633.201479 11864.677062 0.082922 0.075451\n", "\t180 11615.548095 11802.638064 0.079616 0.073662\n", "\t185 11566.309428 11793.285395 0.077288 0.070677\n", "\t190 11575.749336 11763.174836 0.073956 0.068157\n", "\t195 11543.883213 11739.187171 0.071674 0.065917\n", "\t200 11495.598254 11661.501839 0.069762 0.064739\n", "\t205 11490.479524 11650.971928 0.067104 0.062392\n", "\t210 11428.432630 11591.428842 0.065627 0.060987\n", "\t215 11408.594202 11546.502844 0.063437 0.059336\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 11/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.26262e+43 erg / s Luminosity absorbed = 7.78847e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13015.569839 13032.401779 0.682918 0.684003\n", "\t5 12904.875006 12990.461599 0.607993 0.593997\n", "\t10 12842.565429 12898.821890 0.550057 0.540517\n", "\t15 12753.343916 12849.193235 0.505693 0.490550\n", "\t20 12730.133218 12828.953617 0.462479 0.448185\n", "\t25 12664.083226 12759.303136 0.427473 0.414128\n", "\t30 12588.842740 12696.103173 0.399502 0.385663\n", "\t35 12576.400081 12697.730273 0.365150 0.352612\n", "\t40 12571.096532 12662.177207 0.335182 0.324425\n", "\t45 12567.267055 12656.118492 0.308863 0.298253\n", "\t50 12490.762479 12628.623376 0.289858 0.274518\n", "\t55 12473.100822 12623.212204 0.268680 0.252985\n", "\t60 12403.844620 12591.031591 0.252799 0.236744\n", "\t65 12379.017248 12568.788110 0.235858 0.220900\n", "\t70 12384.414417 12504.218648 0.219277 0.209178\n", "\t75 12359.661933 12472.497115 0.205534 0.196505\n", "\t80 12301.567714 12476.433585 0.194798 0.182159\n", "\t85 12282.445145 12444.116662 0.182032 0.170683\n", "\t90 12245.263730 12408.486842 0.172696 0.162078\n", "\t95 12240.211097 12333.668316 0.161126 0.154184\n", "\t100 12166.777597 12303.324181 0.154560 0.145725\n", "\t105 12230.213020 12324.197495 0.141698 0.135636\n", "\t110 12176.872696 12297.477032 0.135522 0.128069\n", "\t115 12173.728933 12298.340501 0.128104 0.120772\n", "\t120 12135.497772 12244.175855 0.121988 0.116145\n", "\t125 12064.394959 12220.660279 0.117167 0.110258\n", "\t130 12010.587301 12203.437058 0.112690 0.104950\n", "\t135 11986.931983 12204.633362 0.107698 0.099736\n", "\t140 12025.244790 12194.368362 0.100791 0.094746\n", "\t145 11995.910253 12188.022851 0.096670 0.090104\n", "\t150 11947.593265 12106.440886 0.093294 0.087625\n", "\t155 11957.542102 12108.357437 0.088421 0.083173\n", "\t160 11910.670130 12056.905371 0.085384 0.080364\n", "\t165 11888.528792 12075.433483 0.081976 0.076282\n", "\t170 11879.429183 12050.882022 0.078285 0.073261\n", "\t175 11864.677062 12025.339383 0.075451 0.070721\n", "\t180 11802.638064 11946.287508 0.073662 0.069445\n", "\t185 11793.285395 11931.799902 0.070677 0.066790\n", "\t190 11763.174836 11906.870608 0.068157 0.064230\n", "\t195 11739.187171 11849.898176 0.065917 0.062563\n", "\t200 11661.501839 11805.166140 0.064739 0.060991\n", "\t205 11650.971928 11750.970412 0.062392 0.059559\n", "\t210 11591.428842 11729.652424 0.060987 0.057636\n", "\t215 11546.502844 11651.688056 0.059336 0.056416\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 12/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.25945e+43 erg / s Luminosity absorbed = 7.81529e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13032.401779 13060.756844 0.684003 0.676972\n", "\t5 12990.461599 13053.655799 0.593997 0.583799\n", "\t10 12898.821890 12940.461814 0.540517 0.534608\n", "\t15 12849.193235 12911.505202 0.490550 0.483174\n", "\t20 12828.953617 12900.206203 0.448185 0.439200\n", "\t25 12759.303136 12863.788022 0.414128 0.400384\n", "\t30 12696.103173 12810.013851 0.385663 0.371252\n", "\t35 12697.730273 12825.414714 0.352612 0.338216\n", "\t40 12662.177207 12801.179422 0.324425 0.310208\n", "\t45 12656.118492 12744.748623 0.298253 0.290688\n", "\t50 12628.623376 12756.940903 0.274518 0.265593\n", "\t55 12623.212204 12719.938507 0.252985 0.246729\n", "\t60 12591.031591 12728.329492 0.236744 0.226363\n", "\t65 12568.788110 12688.708886 0.220900 0.211208\n", "\t70 12504.218648 12660.204209 0.209178 0.197615\n", "\t75 12472.497115 12599.352691 0.196505 0.186919\n", "\t80 12476.433585 12585.793885 0.182159 0.174973\n", "\t85 12444.116662 12589.414222 0.170683 0.162163\n", "\t90 12408.486842 12515.279884 0.162078 0.155383\n", "\t95 12333.668316 12486.897379 0.154184 0.145659\n", "\t100 12303.324181 12430.500086 0.145725 0.138556\n", "\t105 12324.197495 12450.715847 0.135636 0.129425\n", "\t110 12297.477032 12426.647529 0.128069 0.122309\n", "\t115 12298.340501 12435.989959 0.120772 0.114494\n", "\t120 12244.175855 12407.884661 0.116145 0.109221\n", "\t125 12220.660279 12380.981505 0.110258 0.103989\n", "\t130 12203.437058 12351.668125 0.104950 0.099268\n", "\t135 12204.633362 12304.519778 0.099736 0.095737\n", "\t140 12194.368362 12324.560098 0.094746 0.090400\n", "\t145 12188.022851 12265.886677 0.090104 0.087024\n", "\t150 12106.440886 12249.589898 0.087625 0.083028\n", "\t155 12108.357437 12237.977955 0.083173 0.079321\n", "\t160 12056.905371 12207.745803 0.080364 0.076236\n", "\t165 12075.433483 12201.219362 0.076282 0.072831\n", "\t170 12050.882022 12167.349382 0.073261 0.070259\n", "\t175 12025.339383 12095.417592 0.070721 0.068696\n", "\t180 11946.287508 12054.818787 0.069445 0.066319\n", "\t185 11931.799902 12019.446420 0.066790 0.064144\n", "\t190 11906.870608 11990.990568 0.064230 0.061757\n", "\t195 11849.898176 11961.647695 0.062563 0.059858\n", "\t200 11805.166140 11924.727184 0.060991 0.058232\n", "\t205 11750.970412 11858.754302 0.059559 0.056888\n", "\t210 11729.652424 11829.679705 0.057636 0.055114\n", "\t215 11651.688056 11767.828592 0.056416 0.053853\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 13/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.24878e+43 erg / s Luminosity absorbed = 7.90910e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13060.756844 13083.224882 0.676972 0.673697\n", "\t5 13053.655799 13088.802662 0.583799 0.577711\n", "\t10 12940.461814 13027.101067 0.534608 0.522890\n", "\t15 12911.505202 13006.462644 0.483174 0.473946\n", "\t20 12900.206203 12983.888230 0.439200 0.429670\n", "\t25 12863.788022 12971.745863 0.400384 0.388831\n", "\t30 12810.013851 12913.558012 0.371252 0.361218\n", "\t35 12825.414714 12948.249963 0.338216 0.325705\n", "\t40 12801.179422 12912.855169 0.310208 0.301236\n", "\t45 12744.748623 12866.250179 0.290688 0.279866\n", "\t50 12756.940903 12848.652789 0.265593 0.258571\n", "\t55 12719.938507 12848.859656 0.246729 0.236994\n", "\t60 12728.329492 12851.601535 0.226363 0.217969\n", "\t65 12688.708886 12808.658408 0.211208 0.204525\n", "\t70 12660.204209 12799.122574 0.197615 0.189788\n", "\t75 12599.352691 12726.546201 0.186919 0.179897\n", "\t80 12585.793885 12747.782242 0.174973 0.165916\n", "\t85 12589.414222 12750.918368 0.162163 0.154663\n", "\t90 12515.279884 12668.223097 0.155383 0.148464\n", "\t95 12486.897379 12675.791676 0.145659 0.137905\n", "\t100 12430.500086 12620.707010 0.138556 0.130616\n", "\t105 12450.715847 12593.289701 0.129425 0.123022\n", "\t110 12426.647529 12585.676771 0.122309 0.116424\n", "\t115 12435.989959 12537.431630 0.114494 0.110995\n", "\t120 12407.884661 12520.786382 0.109221 0.105135\n", "\t125 12380.981505 12516.716438 0.103989 0.099087\n", "\t130 12351.668125 12472.741826 0.099268 0.095006\n", "\t135 12304.519778 12422.139907 0.095737 0.091867\n", "\t140 12324.560098 12453.435118 0.090400 0.086369\n", "\t145 12265.886677 12389.456714 0.087024 0.083330\n", "\t150 12249.589898 12363.390556 0.083028 0.079554\n", "\t155 12237.977955 12341.590207 0.079321 0.076214\n", "\t160 12207.745803 12291.756528 0.076236 0.073703\n", "\t165 12201.219362 12252.749433 0.072831 0.070982\n", "\t170 12167.349382 12256.286840 0.070259 0.067874\n", "\t175 12095.417592 12182.336237 0.068696 0.066291\n", "\t180 12054.818787 12145.611052 0.066319 0.063949\n", "\t185 12019.446420 12120.021352 0.064144 0.061681\n", "\t190 11990.990568 12095.073710 0.061757 0.059340\n", "\t195 11961.647695 12049.179195 0.059858 0.057897\n", "\t200 11924.727184 11992.337972 0.058232 0.056356\n", "\t205 11858.754302 11949.219275 0.056888 0.054787\n", "\t210 11829.679705 11885.344010 0.055114 0.053671\n", "\t215 11767.828592 11823.561736 0.053853 0.052411\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 14/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.24323e+43 erg / s Luminosity absorbed = 7.95660e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13083.224882 13110.758198 0.673697 0.670180\n", "\t5 13088.802662 13113.028903 0.577711 0.574875\n", "\t10 13027.101067 13093.088462 0.522890 0.515120\n", "\t15 13006.462644 13078.630835 0.473946 0.464486\n", "\t20 12983.888230 13073.991303 0.429670 0.419770\n", "\t25 12971.745863 13110.381168 0.388831 0.375311\n", "\t30 12913.558012 13055.574882 0.361218 0.344822\n", "\t35 12948.249963 13041.575051 0.325705 0.318035\n", "\t40 12912.855169 13022.005042 0.301236 0.292756\n", "\t45 12866.250179 13003.812183 0.279866 0.269639\n", "\t50 12848.652789 12963.032879 0.258571 0.250234\n", "\t55 12848.859656 12982.389906 0.236994 0.229901\n", "\t60 12851.601535 12975.350494 0.217969 0.210624\n", "\t65 12808.658408 12965.368804 0.204525 0.196254\n", "\t70 12799.122574 12890.584285 0.189788 0.185073\n", "\t75 12726.546201 12893.430696 0.179897 0.171994\n", "\t80 12747.782242 12901.177598 0.165916 0.159412\n", "\t85 12750.918368 12894.708333 0.154663 0.148871\n", "\t90 12668.223097 12834.988058 0.148464 0.140906\n", "\t95 12675.791676 12833.456390 0.137905 0.131452\n", "\t100 12620.707010 12807.914014 0.130616 0.123631\n", "\t105 12593.289701 12725.023233 0.123022 0.118612\n", "\t110 12585.676771 12715.728557 0.116424 0.112363\n", "\t115 12537.431630 12681.404241 0.110995 0.106558\n", "\t120 12520.786382 12659.267376 0.105135 0.100885\n", "\t125 12516.716438 12647.301366 0.099087 0.095555\n", "\t130 12472.741826 12601.531043 0.095006 0.091565\n", "\t135 12422.139907 12567.359999 0.091867 0.087506\n", "\t140 12453.435118 12518.962460 0.086369 0.084548\n", "\t145 12389.456714 12478.839839 0.083330 0.080795\n", "\t150 12363.390556 12494.070660 0.079554 0.076008\n", "\t155 12341.590207 12463.099230 0.076214 0.073105\n", "\t160 12291.756528 12383.175568 0.073703 0.071259\n", "\t165 12252.749433 12349.074839 0.070982 0.068442\n", "\t170 12256.286840 12352.713395 0.067874 0.065349\n", "\t175 12182.336237 12260.653843 0.066291 0.064265\n", "\t180 12145.611052 12228.542720 0.063949 0.062020\n", "\t185 12120.021352 12199.166752 0.061681 0.059829\n", "\t190 12095.073710 12177.586624 0.059340 0.057465\n", "\t195 12049.179195 12107.631819 0.057897 0.056255\n", "\t200 11992.337972 12076.794692 0.056356 0.054539\n", "\t205 11949.219275 12033.176341 0.054787 0.053045\n", "\t210 11885.344010 11961.453925 0.053671 0.052030\n", "\t215 11823.561736 11886.901099 0.052411 0.051071\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 15/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.24191e+43 erg / s Luminosity absorbed = 7.96929e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13110.758198 13157.348776 0.670180 0.660872\n", "\t5 13113.028903 13154.537628 0.574875 0.567941\n", "\t10 13093.088462 13122.466176 0.515120 0.508848\n", "\t15 13078.630835 13143.615683 0.464486 0.455821\n", "\t20 13073.991303 13184.013843 0.419770 0.406711\n", "\t25 13110.381168 13204.092069 0.375311 0.367136\n", "\t30 13055.574882 13126.833694 0.344822 0.339156\n", "\t35 13041.575051 13129.174789 0.318035 0.310268\n", "\t40 13022.005042 13109.493361 0.292756 0.285613\n", "\t45 13003.812183 13130.769795 0.269639 0.259172\n", "\t50 12963.032879 13078.921737 0.250234 0.241394\n", "\t55 12982.389906 13088.297094 0.229901 0.222381\n", "\t60 12975.350494 13098.649549 0.210624 0.203764\n", "\t65 12965.368804 13049.630202 0.196254 0.189580\n", "\t70 12890.584285 12990.355074 0.185073 0.179681\n", "\t75 12893.430696 12968.256977 0.171994 0.167990\n", "\t80 12901.177598 12970.293312 0.159412 0.156160\n", "\t85 12894.708333 12968.271549 0.148871 0.145543\n", "\t90 12834.988058 12945.709957 0.140906 0.136494\n", "\t95 12833.456390 12931.106463 0.131452 0.127930\n", "\t100 12807.914014 12919.726451 0.123631 0.120351\n", "\t105 12725.023233 12876.617692 0.118612 0.113582\n", "\t110 12715.728557 12815.972255 0.112363 0.108565\n", "\t115 12681.404241 12832.576329 0.106558 0.102167\n", "\t120 12659.267376 12757.169548 0.100885 0.097927\n", "\t125 12647.301366 12772.943537 0.095555 0.091839\n", "\t130 12601.531043 12694.520268 0.091565 0.088882\n", "\t135 12567.359999 12652.718277 0.087506 0.084953\n", "\t140 12518.962460 12658.308413 0.084548 0.080415\n", "\t145 12478.839839 12630.562968 0.080795 0.077078\n", "\t150 12494.070660 12597.378869 0.076008 0.073622\n", "\t155 12463.099230 12559.021645 0.073105 0.070842\n", "\t160 12383.175568 12498.714833 0.071259 0.068965\n", "\t165 12349.074839 12440.298871 0.068442 0.066375\n", "\t170 12352.713395 12430.528073 0.065349 0.063553\n", "\t175 12260.653843 12372.900884 0.064265 0.062067\n", "\t180 12228.542720 12348.207758 0.062020 0.059442\n", "\t185 12199.166752 12294.218663 0.059829 0.057943\n", "\t190 12177.586624 12281.713874 0.057465 0.055475\n", "\t195 12107.631819 12194.852959 0.056255 0.054501\n", "\t200 12076.794692 12140.252777 0.054539 0.053392\n", "\t205 12033.176341 12102.119966 0.053045 0.051561\n", "\t210 11961.453925 12043.807386 0.052030 0.050461\n", "\t215 11886.901099 11969.967457 0.051071 0.049585\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 16/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.23023e+43 erg / s Luminosity absorbed = 8.07641e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13157.348776 13208.359288 0.660872 0.651688\n", "\t5 13154.537628 13167.132794 0.567941 0.565322\n", "\t10 13122.466176 13185.668397 0.508848 0.502269\n", "\t15 13143.615683 13211.047516 0.455821 0.448037\n", "\t20 13184.013843 13236.736583 0.406711 0.401473\n", "\t25 13204.092069 13206.804495 0.367136 0.366627\n", "\t30 13126.833694 13210.054340 0.339156 0.333981\n", "\t35 13129.174789 13223.266295 0.310268 0.303712\n", "\t40 13109.493361 13144.587177 0.285613 0.283888\n", "\t45 13130.769795 13206.257030 0.259172 0.253767\n", "\t50 13078.921737 13138.179491 0.241394 0.237449\n", "\t55 13088.297094 13162.707185 0.222381 0.217089\n", "\t60 13098.649549 13139.213259 0.203764 0.202532\n", "\t65 13049.630202 13075.253357 0.189580 0.188542\n", "\t70 12990.355074 12998.806255 0.179681 0.178885\n", "\t75 12968.256977 13007.626047 0.167990 0.165915\n", "\t80 12970.293312 12989.924028 0.156160 0.154532\n", "\t85 12968.271549 12950.778512 0.145543 0.144925\n", "\t90 12945.709957 12971.021631 0.136494 0.135208\n", "\t95 12931.106463 12931.545423 0.127930 0.126700\n", "\t100 12919.726451 12967.279881 0.120351 0.118255\n", "\t105 12876.617692 12910.093006 0.113582 0.112558\n", "\t110 12815.972255 12836.599003 0.108565 0.107565\n", "\t115 12832.576329 12807.968042 0.102167 0.102752\n", "\t120 12757.169548 12793.413777 0.097927 0.096426\n", "\t125 12772.943537 12790.716105 0.091839 0.091588\n", "\t130 12694.520268 12762.147943 0.088882 0.086533\n", "\t135 12652.718277 12707.570217 0.084953 0.083250\n", "\t140 12658.308413 12723.808197 0.080415 0.078572\n", "\t145 12630.562968 12631.971267 0.077078 0.076837\n", "\t150 12597.378869 12613.669514 0.073622 0.072776\n", "\t155 12559.021645 12571.136723 0.070842 0.070026\n", "\t160 12498.714833 12582.788478 0.068965 0.066758\n", "\t165 12440.298871 12515.372784 0.066375 0.064675\n", "\t170 12430.528073 12463.372449 0.063553 0.062916\n", "\t175 12372.900884 12414.922501 0.062067 0.060897\n", "\t180 12348.207758 12374.558711 0.059442 0.058655\n", "\t185 12294.218663 12325.073521 0.057943 0.056945\n", "\t190 12281.713874 12271.369723 0.055475 0.055274\n", "\t195 12194.852959 12215.455284 0.054501 0.054039\n", "\t200 12140.252777 12184.530283 0.053392 0.052251\n", "\t205 12102.119966 12150.063082 0.051561 0.050542\n", "\t210 12043.807386 12106.155134 0.050461 0.049397\n", "\t215 11969.967457 12017.087711 0.049585 0.048618\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 17/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.23593e+43 erg / s Luminosity absorbed = 8.01895e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13208.359288 13180.535471 0.651688 0.658012\n", "\t5 13167.132794 13197.522520 0.565322 0.562610\n", "\t10 13185.668397 13174.694866 0.502269 0.505905\n", "\t15 13211.047516 13218.240350 0.448037 0.447772\n", "\t20 13236.736583 13241.225892 0.401473 0.403604\n", "\t25 13206.804495 13281.145141 0.366627 0.360381\n", "\t30 13210.054340 13267.313536 0.333981 0.328384\n", "\t35 13223.266295 13286.738073 0.303712 0.297458\n", "\t40 13144.587177 13181.026751 0.283888 0.279818\n", "\t45 13206.257030 13250.028942 0.253767 0.251366\n", "\t50 13138.179491 13221.355944 0.237449 0.232221\n", "\t55 13162.707185 13184.596993 0.217089 0.216556\n", "\t60 13139.213259 13167.869456 0.202532 0.200682\n", "\t65 13075.253357 13157.266030 0.188542 0.184969\n", "\t70 12998.806255 13067.082659 0.178885 0.176158\n", "\t75 13007.626047 13063.765451 0.165915 0.163820\n", "\t80 12989.924028 13107.403255 0.154532 0.150006\n", "\t85 12950.778512 13019.626254 0.144925 0.142593\n", "\t90 12971.021631 13075.761908 0.135208 0.131608\n", "\t95 12931.545423 12974.276546 0.126700 0.125865\n", "\t100 12967.279881 12999.965074 0.118255 0.117449\n", "\t105 12910.093006 12960.919781 0.112558 0.111632\n", "\t110 12836.599003 12862.321828 0.107565 0.106997\n", "\t115 12807.968042 12888.724834 0.102752 0.100015\n", "\t120 12793.413777 12843.322883 0.096426 0.094920\n", "\t125 12790.716105 12840.666968 0.091588 0.089666\n", "\t130 12762.147943 12747.676705 0.086533 0.087161\n", "\t135 12707.570217 12724.372823 0.083250 0.083115\n", "\t140 12723.808197 12748.930414 0.078572 0.078337\n", "\t145 12631.971267 12717.078441 0.076837 0.074678\n", "\t150 12613.669514 12663.247909 0.072776 0.071613\n", "\t155 12571.136723 12639.626393 0.070026 0.069149\n", "\t160 12582.788478 12605.013851 0.066758 0.066224\n", "\t165 12515.372784 12529.037411 0.064675 0.064437\n", "\t170 12463.372449 12512.441055 0.062916 0.061855\n", "\t175 12414.922501 12504.982434 0.060897 0.059360\n", "\t180 12374.558711 12421.316875 0.058655 0.057855\n", "\t185 12325.073521 12407.198976 0.056945 0.055458\n", "\t190 12271.369723 12340.725468 0.055274 0.054148\n", "\t195 12215.455284 12284.031240 0.054039 0.052758\n", "\t200 12184.530283 12240.133374 0.052251 0.051248\n", "\t205 12150.063082 12197.305230 0.050542 0.049849\n", "\t210 12106.155134 12155.057432 0.049397 0.048580\n", "\t215 12017.087711 12048.113577 0.048618 0.048016\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 18/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.23486e+43 erg / s Luminosity absorbed = 8.02956e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13180.535471 13166.826599 0.658012 0.661773\n", "\t5 13197.522520 13161.225439 0.562610 0.570444\n", "\t10 13174.694866 13164.033321 0.505905 0.508207\n", "\t15 13218.240350 13201.849346 0.447772 0.450656\n", "\t20 13241.225892 13207.232370 0.403604 0.408133\n", "\t25 13281.145141 13265.773013 0.360381 0.360941\n", "\t30 13267.313536 13256.622272 0.328384 0.330341\n", "\t35 13286.738073 13303.315390 0.297458 0.296541\n", "\t40 13181.026751 13239.259991 0.279818 0.276023\n", "\t45 13250.028942 13284.146284 0.251366 0.251606\n", "\t50 13221.355944 13281.863495 0.232221 0.229871\n", "\t55 13184.596993 13291.777175 0.216556 0.210540\n", "\t60 13167.869456 13273.989481 0.200682 0.195987\n", "\t65 13157.266030 13234.724167 0.184969 0.181474\n", "\t70 13067.082659 13073.022860 0.176158 0.176104\n", "\t75 13063.765451 13104.206567 0.163820 0.161397\n", "\t80 13107.403255 13144.009619 0.150006 0.148310\n", "\t85 13019.626254 13095.391976 0.142593 0.139576\n", "\t90 13075.761908 13127.347738 0.131608 0.129448\n", "\t95 12974.276546 13047.484719 0.125865 0.123275\n", "\t100 12999.965074 13054.715390 0.117449 0.115835\n", "\t105 12960.919781 13030.075993 0.111632 0.109099\n", "\t110 12862.321828 12933.735701 0.106997 0.104811\n", "\t115 12888.724834 12932.599283 0.100015 0.099248\n", "\t120 12843.322883 12887.477861 0.094920 0.094004\n", "\t125 12840.666968 12856.250561 0.089666 0.089628\n", "\t130 12747.676705 12868.962924 0.087161 0.083908\n", "\t135 12724.372823 12788.877268 0.083115 0.081364\n", "\t140 12748.930414 12795.467458 0.078337 0.077004\n", "\t145 12717.078441 12763.950779 0.074678 0.073739\n", "\t150 12663.247909 12718.648984 0.071613 0.070597\n", "\t155 12639.626393 12715.537044 0.069149 0.067710\n", "\t160 12605.013851 12650.932274 0.066224 0.065329\n", "\t165 12529.037411 12607.924877 0.064437 0.062762\n", "\t170 12512.441055 12608.704017 0.061855 0.059987\n", "\t175 12504.982434 12586.366176 0.059360 0.057683\n", "\t180 12421.316875 12528.864095 0.057855 0.055953\n", "\t185 12407.198976 12477.716927 0.055458 0.054166\n", "\t190 12340.725468 12414.012399 0.054148 0.052969\n", "\t195 12284.031240 12346.489093 0.052758 0.051603\n", "\t200 12240.133374 12285.180423 0.051248 0.050515\n", "\t205 12197.305230 12227.179556 0.049849 0.049302\n", "\t210 12155.057432 12166.263653 0.048580 0.048302\n", "\t215 12048.113577 12071.857550 0.048016 0.047522\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 19/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22794e+43 erg / s Luminosity absorbed = 8.08940e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13166.826599 13195.519766 0.661773 0.655923\n", "\t5 13161.225439 13210.741764 0.570444 0.561298\n", "\t10 13164.033321 13235.889833 0.508207 0.496477\n", "\t15 13201.849346 13257.976739 0.450656 0.442108\n", "\t20 13207.232370 13244.393460 0.408133 0.404152\n", "\t25 13265.773013 13296.492302 0.360941 0.360900\n", "\t30 13256.622272 13297.435443 0.330341 0.328098\n", "\t35 13303.315390 13345.323924 0.296541 0.294771\n", "\t40 13239.259991 13291.955653 0.276023 0.273309\n", "\t45 13284.146284 13330.280434 0.251606 0.248861\n", "\t50 13281.863495 13308.757907 0.229871 0.230862\n", "\t55 13291.777175 13304.285435 0.210540 0.211845\n", "\t60 13273.989481 13296.137731 0.195987 0.196026\n", "\t65 13234.724167 13270.557153 0.181474 0.180536\n", "\t70 13073.022860 13177.400597 0.176104 0.171785\n", "\t75 13104.206567 13183.943549 0.161397 0.158635\n", "\t80 13144.009619 13178.261898 0.148310 0.148387\n", "\t85 13095.391976 13152.438819 0.139576 0.138476\n", "\t90 13127.347738 13215.251923 0.129448 0.126976\n", "\t95 13047.484719 13142.709904 0.123275 0.121220\n", "\t100 13054.715390 13117.686046 0.115835 0.114414\n", "\t105 13030.075993 13084.549148 0.109099 0.107408\n", "\t110 12933.735701 13016.108830 0.104811 0.103248\n", "\t115 12932.599283 13034.016323 0.099248 0.095773\n", "\t120 12887.477861 12996.130306 0.094004 0.091197\n", "\t125 12856.250561 12941.186178 0.089628 0.087218\n", "\t130 12868.962924 12925.433879 0.083908 0.082603\n", "\t135 12788.877268 12880.728193 0.081364 0.079023\n", "\t140 12795.467458 12847.699135 0.077004 0.075342\n", "\t145 12763.950779 12794.111032 0.073739 0.072687\n", "\t150 12718.648984 12736.936646 0.070597 0.069979\n", "\t155 12715.537044 12710.772364 0.067710 0.067228\n", "\t160 12650.932274 12649.594065 0.065329 0.065043\n", "\t165 12607.924877 12622.622036 0.062762 0.062659\n", "\t170 12608.704017 12611.817207 0.059987 0.059724\n", "\t175 12586.366176 12585.348534 0.057683 0.057595\n", "\t180 12528.864095 12550.454375 0.055953 0.055510\n", "\t185 12477.716927 12478.754739 0.054166 0.053976\n", "\t190 12414.012399 12444.539262 0.052969 0.052245\n", "\t195 12346.489093 12373.740602 0.051603 0.050962\n", "\t200 12285.180423 12313.701285 0.050515 0.049849\n", "\t205 12227.179556 12241.245001 0.049302 0.048828\n", "\t210 12166.263653 12183.394427 0.048302 0.047714\n", "\t215 12071.857550 12099.018816 0.047522 0.046976\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 20/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22565e+43 erg / s Luminosity absorbed = 8.10622e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13195.519766 13218.392110 0.655923 0.649349\n", "\t5 13210.741764 13203.590883 0.561298 0.562997\n", "\t10 13235.889833 13261.694418 0.496477 0.494806\n", "\t15 13257.976739 13298.064876 0.442108 0.438235\n", "\t20 13244.393460 13250.974631 0.404152 0.404452\n", "\t25 13296.492302 13305.202082 0.360900 0.360661\n", "\t30 13297.435443 13329.187270 0.328098 0.326259\n", "\t35 13345.323924 13351.564942 0.294771 0.295542\n", "\t40 13291.955653 13350.505065 0.273309 0.269488\n", "\t45 13330.280434 13381.278649 0.248861 0.244352\n", "\t50 13308.757907 13350.401391 0.230862 0.227866\n", "\t55 13304.285435 13322.777854 0.211845 0.212514\n", "\t60 13296.137731 13338.479356 0.196026 0.194321\n", "\t65 13270.557153 13301.899148 0.180536 0.180111\n", "\t70 13177.400597 13281.385654 0.171785 0.167801\n", "\t75 13183.943549 13294.307652 0.158635 0.154838\n", "\t80 13178.261898 13267.966077 0.148387 0.145283\n", "\t85 13152.438819 13231.259161 0.138476 0.135669\n", "\t90 13215.251923 13243.603339 0.126976 0.126946\n", "\t95 13142.709904 13245.685361 0.121220 0.117504\n", "\t100 13117.686046 13177.508230 0.114414 0.112469\n", "\t105 13084.549148 13166.889212 0.107408 0.105195\n", "\t110 13016.108830 13166.758929 0.103248 0.099258\n", "\t115 13034.016323 13147.285482 0.095773 0.092737\n", "\t120 12996.130306 13132.027948 0.091197 0.087745\n", "\t125 12941.186178 13024.016775 0.087218 0.085256\n", "\t130 12925.433879 13018.325987 0.082603 0.080448\n", "\t135 12880.728193 13030.706325 0.079023 0.075968\n", "\t140 12847.699135 12937.788695 0.075342 0.073766\n", "\t145 12794.111032 12911.995175 0.072687 0.069983\n", "\t150 12736.936646 12823.855253 0.069979 0.068040\n", "\t155 12710.772364 12801.548304 0.067228 0.065332\n", "\t160 12649.594065 12716.134743 0.065043 0.063471\n", "\t165 12622.622036 12684.472936 0.062659 0.061499\n", "\t170 12611.817207 12681.171176 0.059724 0.058254\n", "\t175 12585.348534 12622.297246 0.057595 0.056786\n", "\t180 12550.454375 12571.964450 0.055510 0.054959\n", "\t185 12478.754739 12526.205063 0.053976 0.053273\n", "\t190 12444.539262 12472.393760 0.052245 0.051593\n", "\t195 12373.740602 12392.121290 0.050962 0.050498\n", "\t200 12313.701285 12324.192334 0.049849 0.049509\n", "\t205 12241.245001 12255.227451 0.048828 0.048525\n", "\t210 12183.394427 12234.356154 0.047714 0.046820\n", "\t215 12099.018816 12116.932757 0.046976 0.046558\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 21/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22962e+43 erg / s Luminosity absorbed = 8.07205e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13218.392110 13193.465240 0.649349 0.653651\n", "\t5 13203.590883 13227.367950 0.562997 0.559593\n", "\t10 13261.694418 13274.159820 0.494806 0.492442\n", "\t15 13298.064876 13293.745040 0.438235 0.439388\n", "\t20 13250.974631 13288.270014 0.404452 0.399678\n", "\t25 13305.202082 13292.192089 0.360661 0.361950\n", "\t30 13329.187270 13317.502835 0.326259 0.325734\n", "\t35 13351.564942 13401.993759 0.295542 0.290505\n", "\t40 13350.505065 13359.932654 0.269488 0.270203\n", "\t45 13381.278649 13407.014769 0.244352 0.243462\n", "\t50 13350.401391 13409.046073 0.227866 0.224079\n", "\t55 13322.777854 13340.715423 0.212514 0.210605\n", "\t60 13338.479356 13381.415152 0.194321 0.191791\n", "\t65 13301.899148 13377.188342 0.180111 0.176521\n", "\t70 13281.385654 13334.189209 0.167801 0.165216\n", "\t75 13294.307652 13348.055264 0.154838 0.152572\n", "\t80 13267.966077 13322.318785 0.145283 0.142626\n", "\t85 13231.259161 13272.891041 0.135669 0.133746\n", "\t90 13243.603339 13280.211749 0.126946 0.125982\n", "\t95 13245.685361 13267.189239 0.117504 0.116621\n", "\t100 13177.508230 13215.630844 0.112469 0.111315\n", "\t105 13166.889212 13175.010755 0.105195 0.105045\n", "\t110 13166.758929 13200.228786 0.099258 0.097975\n", "\t115 13147.285482 13162.211276 0.092737 0.092293\n", "\t120 13132.027948 13106.849919 0.087745 0.088225\n", "\t125 13024.016775 13055.657025 0.085256 0.084607\n", "\t130 13018.325987 13012.306316 0.080448 0.080817\n", "\t135 13030.706325 12997.425810 0.075968 0.076277\n", "\t140 12937.788695 12945.009959 0.073766 0.073460\n", "\t145 12911.995175 12918.379242 0.069983 0.070102\n", "\t150 12823.855253 12811.234197 0.068040 0.068201\n", "\t155 12801.548304 12773.399696 0.065332 0.066079\n", "\t160 12716.134743 12751.645015 0.063471 0.063025\n", "\t165 12684.472936 12724.598687 0.061499 0.060578\n", "\t170 12681.171176 12676.535769 0.058254 0.058497\n", "\t175 12622.297246 12617.758328 0.056786 0.056922\n", "\t180 12571.964450 12566.081512 0.054959 0.055095\n", "\t185 12526.205063 12519.386548 0.053273 0.053572\n", "\t190 12472.393760 12501.909282 0.051593 0.051469\n", "\t195 12392.121290 12437.036377 0.050498 0.049955\n", "\t200 12324.192334 12352.874273 0.049509 0.049206\n", "\t205 12255.227451 12282.392519 0.048525 0.048163\n", "\t210 12234.356154 12217.553356 0.046820 0.047145\n", "\t215 12116.932757 12106.194987 0.046558 0.046731\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 22/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22206e+43 erg / s Luminosity absorbed = 8.13794e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13193.465240 13195.466506 0.653651 0.655681\n", "\t5 13227.367950 13234.149938 0.559593 0.558807\n", "\t10 13274.159820 13250.119206 0.492442 0.497110\n", "\t15 13293.745040 13283.570386 0.439388 0.443546\n", "\t20 13288.270014 13285.196098 0.399678 0.401518\n", "\t25 13292.192089 13297.903037 0.361950 0.363788\n", "\t30 13317.502835 13304.903433 0.325734 0.330382\n", "\t35 13401.993759 13348.630577 0.290505 0.296154\n", "\t40 13359.932654 13373.208881 0.270203 0.268559\n", "\t45 13407.014769 13409.696334 0.243462 0.243838\n", "\t50 13409.046073 13446.583648 0.224079 0.221260\n", "\t55 13340.715423 13363.747733 0.210605 0.209475\n", "\t60 13381.415152 13445.718362 0.191791 0.188191\n", "\t65 13377.188342 13451.866500 0.176521 0.173362\n", "\t70 13334.189209 13420.709807 0.165216 0.161404\n", "\t75 13348.055264 13404.576966 0.152572 0.150490\n", "\t80 13322.318785 13376.205532 0.142626 0.141152\n", "\t85 13272.891041 13367.313511 0.133746 0.130660\n", "\t90 13280.211749 13362.190573 0.125982 0.122816\n", "\t95 13267.189239 13315.019271 0.116621 0.114927\n", "\t100 13215.630844 13284.482648 0.111315 0.108049\n", "\t105 13175.010755 13228.775249 0.105045 0.103409\n", "\t110 13200.228786 13250.373347 0.097975 0.096366\n", "\t115 13162.211276 13163.949582 0.092293 0.092277\n", "\t120 13106.849919 13195.199414 0.088225 0.086233\n", "\t125 13055.657025 13168.769528 0.084607 0.081787\n", "\t130 13012.306316 13085.969357 0.080817 0.079046\n", "\t135 12997.425810 13075.101627 0.076277 0.074932\n", "\t140 12945.009959 12966.424621 0.073460 0.073250\n", "\t145 12918.379242 12951.468346 0.070102 0.069209\n", "\t150 12811.234197 12890.777445 0.068201 0.066757\n", "\t155 12773.399696 12846.206270 0.066079 0.064555\n", "\t160 12751.645015 12803.339348 0.063025 0.062149\n", "\t165 12724.598687 12797.147056 0.060578 0.059464\n", "\t170 12676.535769 12797.437041 0.058497 0.056845\n", "\t175 12617.758328 12708.051794 0.056922 0.055482\n", "\t180 12566.081512 12637.715416 0.055095 0.054127\n", "\t185 12519.386548 12574.299005 0.053572 0.052499\n", "\t190 12501.909282 12531.863048 0.051469 0.050803\n", "\t195 12437.036377 12496.462195 0.049955 0.049082\n", "\t200 12352.874273 12416.783340 0.049206 0.048305\n", "\t205 12282.392519 12349.739513 0.048163 0.047137\n", "\t210 12217.553356 12254.084604 0.047145 0.046408\n", "\t215 12106.194987 12160.167476 0.046731 0.045780\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 23/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22210e+43 erg / s Luminosity absorbed = 8.13819e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13195.466506 13216.402157 0.655681 0.652541\n", "\t5 13234.149938 13237.225427 0.558807 0.559506\n", "\t10 13250.119206 13252.209818 0.497110 0.497424\n", "\t15 13283.570386 13243.043291 0.443546 0.449203\n", "\t20 13285.196098 13279.986369 0.401518 0.401513\n", "\t25 13297.903037 13307.974117 0.363788 0.360489\n", "\t30 13304.903433 13330.259306 0.330382 0.326376\n", "\t35 13348.630577 13361.962084 0.296154 0.295078\n", "\t40 13373.208881 13394.308987 0.268559 0.266887\n", "\t45 13409.696334 13436.815696 0.243838 0.242688\n", "\t50 13446.583648 13461.659991 0.221260 0.221027\n", "\t55 13363.747733 13443.756797 0.209475 0.205664\n", "\t60 13445.718362 13439.851226 0.188191 0.188685\n", "\t65 13451.866500 13447.748126 0.173362 0.174022\n", "\t70 13420.709807 13423.536597 0.161404 0.161184\n", "\t75 13404.576966 13423.527026 0.150490 0.149976\n", "\t80 13376.205532 13431.673399 0.141152 0.138622\n", "\t85 13367.313511 13367.261282 0.130660 0.130878\n", "\t90 13362.190573 13392.504047 0.122816 0.121643\n", "\t95 13315.019271 13379.838946 0.114927 0.113123\n", "\t100 13284.482648 13305.320906 0.108049 0.107698\n", "\t105 13228.775249 13299.509838 0.103409 0.101494\n", "\t110 13250.373347 13303.155862 0.096366 0.095391\n", "\t115 13163.949582 13239.401607 0.092277 0.090675\n", "\t120 13195.199414 13225.853168 0.086233 0.085758\n", "\t125 13168.769528 13182.779625 0.081787 0.082242\n", "\t130 13085.969357 13158.856980 0.079046 0.077699\n", "\t135 13075.101627 13080.616663 0.074932 0.075268\n", "\t140 12966.424621 13027.168954 0.073250 0.072122\n", "\t145 12951.468346 13000.819102 0.069209 0.068449\n", "\t150 12890.777445 12985.503070 0.066757 0.065123\n", "\t155 12846.206270 12930.672063 0.064555 0.062974\n", "\t160 12803.339348 12854.302745 0.062149 0.061227\n", "\t165 12797.147056 12825.373038 0.059464 0.058748\n", "\t170 12797.437041 12819.415527 0.056845 0.056058\n", "\t175 12708.051794 12708.497312 0.055482 0.055067\n", "\t180 12637.715416 12639.246258 0.054127 0.054019\n", "\t185 12574.299005 12631.362226 0.052499 0.051403\n", "\t190 12531.863048 12559.773216 0.050803 0.050265\n", "\t195 12496.462195 12612.585828 0.049082 0.047577\n", "\t200 12416.783340 12489.714343 0.048305 0.047266\n", "\t205 12349.739513 12414.869027 0.047137 0.046149\n", "\t210 12254.084604 12333.119968 0.046408 0.045325\n", "\t215 12160.167476 12218.394834 0.045780 0.044907\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 24/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21424e+43 erg / s Luminosity absorbed = 8.21243e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13216.402157 13214.412341 0.652541 0.654355\n", "\t5 13237.225427 13266.546018 0.559506 0.557950\n", "\t10 13252.209818 13266.925438 0.497424 0.495958\n", "\t15 13243.043291 13321.265822 0.449203 0.439403\n", "\t20 13279.986369 13300.999426 0.401513 0.399679\n", "\t25 13307.974117 13340.522690 0.360489 0.359001\n", "\t30 13330.259306 13360.096557 0.326376 0.324063\n", "\t35 13361.962084 13425.244338 0.295078 0.290336\n", "\t40 13394.308987 13406.012752 0.266887 0.267246\n", "\t45 13436.815696 13470.738156 0.242688 0.239683\n", "\t50 13461.659991 13514.566779 0.221027 0.216757\n", "\t55 13443.756797 13446.809708 0.205664 0.204399\n", "\t60 13439.851226 13495.051511 0.188685 0.185571\n", "\t65 13447.748126 13478.049586 0.174022 0.171728\n", "\t70 13423.536597 13478.243438 0.161184 0.158297\n", "\t75 13423.527026 13447.494233 0.149976 0.148463\n", "\t80 13431.673399 13436.926051 0.138622 0.138342\n", "\t85 13367.261282 13405.984926 0.130878 0.129554\n", "\t90 13392.504047 13411.477844 0.121643 0.120318\n", "\t95 13379.838946 13397.359911 0.113123 0.112874\n", "\t100 13305.320906 13360.714798 0.107698 0.106382\n", "\t105 13299.509838 13331.840059 0.101494 0.100378\n", "\t110 13303.155862 13283.695143 0.095391 0.095469\n", "\t115 13239.401607 13269.496703 0.090675 0.089693\n", "\t120 13225.853168 13216.630334 0.085758 0.085353\n", "\t125 13182.779625 13188.100904 0.082242 0.081391\n", "\t130 13158.856980 13171.671286 0.077699 0.077299\n", "\t135 13080.616663 13114.333697 0.075268 0.074280\n", "\t140 13027.168954 13087.991186 0.072122 0.070821\n", "\t145 13000.819102 13030.686619 0.068449 0.067675\n", "\t150 12985.503070 12997.630280 0.065123 0.064682\n", "\t155 12930.672063 12953.381884 0.062974 0.062689\n", "\t160 12854.302745 12869.969481 0.061227 0.060758\n", "\t165 12825.373038 12850.813423 0.058748 0.058063\n", "\t170 12819.415527 12822.577348 0.056058 0.055673\n", "\t175 12708.497312 12719.584918 0.055067 0.054788\n", "\t180 12639.246258 12641.476771 0.054019 0.053491\n", "\t185 12631.362226 12625.848643 0.051403 0.051273\n", "\t190 12559.773216 12633.826899 0.050265 0.049295\n", "\t195 12612.585828 12568.263736 0.047577 0.047966\n", "\t200 12489.714343 12497.684613 0.047266 0.046947\n", "\t205 12414.869027 12403.984719 0.046149 0.046160\n", "\t210 12333.119968 12336.380238 0.045325 0.045221\n", "\t215 12218.394834 12223.584261 0.044907 0.044754\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 25/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21658e+43 erg / s Luminosity absorbed = 8.19057e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13214.412341 13226.637262 0.654355 0.649960\n", "\t5 13266.546018 13273.317748 0.557950 0.557229\n", "\t10 13266.925438 13308.672690 0.495958 0.489640\n", "\t15 13321.265822 13362.955586 0.439403 0.433103\n", "\t20 13300.999426 13295.175381 0.399679 0.399624\n", "\t25 13340.522690 13344.623968 0.359001 0.358202\n", "\t30 13360.096557 13391.442684 0.324063 0.322013\n", "\t35 13425.244338 13466.503250 0.290336 0.288319\n", "\t40 13406.012752 13436.991080 0.267246 0.265779\n", "\t45 13470.738156 13513.082253 0.239683 0.237983\n", "\t50 13514.566779 13499.586715 0.216757 0.218713\n", "\t55 13446.809708 13471.751080 0.204399 0.203875\n", "\t60 13495.051511 13474.365905 0.185571 0.187673\n", "\t65 13478.049586 13495.089584 0.171728 0.171978\n", "\t70 13478.243438 13506.466553 0.158297 0.158035\n", "\t75 13447.494233 13472.618979 0.148463 0.147689\n", "\t80 13436.926051 13426.092569 0.138342 0.138483\n", "\t85 13405.984926 13435.279967 0.129554 0.127994\n", "\t90 13411.477844 13398.549637 0.120318 0.120507\n", "\t95 13397.359911 13392.899837 0.112874 0.112153\n", "\t100 13360.714798 13355.280323 0.106382 0.106066\n", "\t105 13331.840059 13350.876967 0.100378 0.100158\n", "\t110 13283.695143 13262.987994 0.095469 0.096234\n", "\t115 13269.496703 13269.261900 0.089693 0.089592\n", "\t120 13216.630334 13220.271865 0.085353 0.085028\n", "\t125 13188.100904 13227.413382 0.081391 0.080367\n", "\t130 13171.671286 13183.623645 0.077299 0.076439\n", "\t135 13114.333697 13154.819805 0.074280 0.072986\n", "\t140 13087.991186 13073.616423 0.070821 0.070857\n", "\t145 13030.686619 13043.634237 0.067675 0.067414\n", "\t150 12997.630280 12994.570180 0.064682 0.064596\n", "\t155 12953.381884 12946.036078 0.062689 0.062456\n", "\t160 12869.969481 12899.812989 0.060758 0.060078\n", "\t165 12850.813423 12860.432622 0.058063 0.057899\n", "\t170 12822.577348 12758.171162 0.055673 0.056694\n", "\t175 12719.584918 12736.508057 0.054788 0.054428\n", "\t180 12641.476771 12693.757431 0.053491 0.052712\n", "\t185 12625.848643 12627.684428 0.051273 0.051319\n", "\t190 12633.826899 12605.649726 0.049295 0.049639\n", "\t195 12568.263736 12519.805078 0.047966 0.048512\n", "\t200 12497.684613 12441.375388 0.046947 0.047643\n", "\t205 12403.984719 12390.758682 0.046160 0.046193\n", "\t210 12336.380238 12341.490806 0.045221 0.045034\n", "\t215 12223.584261 12269.951528 0.044754 0.044087\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 26/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.22024e+43 erg / s Luminosity absorbed = 8.15492e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13226.637262 13224.038802 0.649960 0.649671\n", "\t5 13273.317748 13254.407087 0.557229 0.561481\n", "\t10 13308.672690 13284.603566 0.489640 0.492310\n", "\t15 13362.955586 13302.545105 0.433103 0.439342\n", "\t20 13295.175381 13292.262849 0.399624 0.398956\n", "\t25 13344.623968 13342.035745 0.358202 0.357113\n", "\t30 13391.442684 13395.651625 0.322013 0.320898\n", "\t35 13466.503250 13464.557456 0.288319 0.289676\n", "\t40 13436.991080 13440.515426 0.265779 0.266129\n", "\t45 13513.082253 13515.022982 0.237983 0.237762\n", "\t50 13499.586715 13510.184201 0.218713 0.218257\n", "\t55 13471.751080 13504.086537 0.203875 0.201539\n", "\t60 13474.365905 13468.713278 0.187673 0.187599\n", "\t65 13495.089584 13499.187263 0.171978 0.172129\n", "\t70 13506.466553 13503.375732 0.158035 0.158278\n", "\t75 13472.618979 13471.049950 0.147689 0.147953\n", "\t80 13426.092569 13469.814316 0.138483 0.137510\n", "\t85 13435.279967 13466.753912 0.127994 0.127882\n", "\t90 13398.549637 13414.734938 0.120507 0.121051\n", "\t95 13392.899837 13424.523000 0.112153 0.112245\n", "\t100 13355.280323 13380.318364 0.106066 0.106457\n", "\t105 13350.876967 13352.257978 0.100158 0.100303\n", "\t110 13262.987994 13320.346400 0.096234 0.094864\n", "\t115 13269.261900 13299.997262 0.089592 0.089149\n", "\t120 13220.271865 13232.512535 0.085028 0.084875\n", "\t125 13227.413382 13242.894290 0.080367 0.080056\n", "\t130 13183.623645 13202.446665 0.076439 0.076082\n", "\t135 13154.819805 13156.812152 0.072986 0.073040\n", "\t140 13073.616423 13076.700422 0.070857 0.070527\n", "\t145 13043.634237 13053.722896 0.067414 0.067399\n", "\t150 12994.570180 13006.987684 0.064596 0.064253\n", "\t155 12946.036078 12967.210224 0.062456 0.062075\n", "\t160 12899.812989 12945.463999 0.060078 0.059531\n", "\t165 12860.432622 12854.793486 0.057899 0.057858\n", "\t170 12758.171162 12796.353980 0.056694 0.056173\n", "\t175 12736.508057 12792.868597 0.054428 0.053488\n", "\t180 12693.757431 12696.129047 0.052712 0.052492\n", "\t185 12627.684428 12656.106175 0.051319 0.050890\n", "\t190 12605.649726 12598.256389 0.049639 0.049612\n", "\t195 12519.805078 12518.316319 0.048512 0.048575\n", "\t200 12441.375388 12456.817563 0.047643 0.047361\n", "\t205 12390.758682 12358.693871 0.046193 0.046701\n", "\t210 12341.490806 12305.940437 0.045034 0.045632\n", "\t215 12269.951528 12229.092749 0.044087 0.044648\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 27/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21669e+43 erg / s Luminosity absorbed = 8.18670e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13224.038802 13227.901174 0.649671 0.649570\n", "\t5 13254.407087 13253.437801 0.561481 0.562809\n", "\t10 13284.603566 13315.822899 0.492310 0.487354\n", "\t15 13302.545105 13308.504973 0.439342 0.441817\n", "\t20 13292.262849 13313.086028 0.398956 0.398092\n", "\t25 13342.035745 13362.257074 0.357113 0.355488\n", "\t30 13395.651625 13385.926898 0.320898 0.322538\n", "\t35 13464.557456 13455.152434 0.289676 0.290891\n", "\t40 13440.515426 13468.663938 0.266129 0.263483\n", "\t45 13515.022982 13519.866805 0.237762 0.238079\n", "\t50 13510.184201 13506.898673 0.218257 0.219233\n", "\t55 13504.086537 13549.458452 0.201539 0.199181\n", "\t60 13468.713278 13509.018362 0.187599 0.185527\n", "\t65 13499.187263 13486.072442 0.172129 0.172556\n", "\t70 13503.375732 13525.678046 0.158278 0.157214\n", "\t75 13471.049950 13482.036613 0.147953 0.146788\n", "\t80 13469.814316 13503.408657 0.137510 0.135660\n", "\t85 13466.753912 13518.750962 0.127882 0.125845\n", "\t90 13414.734938 13448.183311 0.121051 0.118639\n", "\t95 13424.523000 13466.878706 0.112245 0.110304\n", "\t100 13380.318364 13393.365097 0.106457 0.105151\n", "\t105 13352.257978 13320.206391 0.100303 0.100640\n", "\t110 13320.346400 13296.148156 0.094864 0.095177\n", "\t115 13299.997262 13277.749438 0.089149 0.089847\n", "\t120 13232.512535 13255.594266 0.084875 0.084687\n", "\t125 13242.894290 13216.769722 0.080056 0.080475\n", "\t130 13202.446665 13134.280941 0.076082 0.077466\n", "\t135 13156.812152 13143.270696 0.073040 0.073214\n", "\t140 13076.700422 13118.062526 0.070527 0.069721\n", "\t145 13053.722896 13129.778037 0.067399 0.066201\n", "\t150 13006.987684 13016.758140 0.064253 0.064372\n", "\t155 12967.210224 12952.608913 0.062075 0.062374\n", "\t160 12945.463999 12972.195659 0.059531 0.058840\n", "\t165 12854.793486 12901.782640 0.057858 0.057205\n", "\t170 12796.353980 12836.110004 0.056173 0.055623\n", "\t175 12792.868597 12813.578279 0.053488 0.053283\n", "\t180 12696.129047 12746.183324 0.052492 0.052010\n", "\t185 12656.106175 12642.400168 0.050890 0.051181\n", "\t190 12598.256389 12623.518838 0.049612 0.049208\n", "\t195 12518.316319 12545.983504 0.048575 0.048173\n", "\t200 12456.817563 12486.360144 0.047361 0.046895\n", "\t205 12358.693871 12395.242697 0.046701 0.046184\n", "\t210 12305.940437 12343.814221 0.045632 0.045038\n", "\t215 12229.092749 12244.749981 0.044648 0.044401\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 28/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21499e+43 erg / s Luminosity absorbed = 8.20205e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13227.901174 13233.150013 0.649570 0.650204\n", "\t5 13253.437801 13246.153944 0.562809 0.561930\n", "\t10 13315.822899 13321.909830 0.487354 0.489863\n", "\t15 13308.504973 13364.961838 0.441817 0.435686\n", "\t20 13313.086028 13391.543417 0.398092 0.390502\n", "\t25 13362.257074 13380.970341 0.355488 0.354792\n", "\t30 13385.926898 13424.798756 0.322538 0.319017\n", "\t35 13455.152434 13461.928210 0.290891 0.291085\n", "\t40 13468.663938 13503.499504 0.263483 0.261981\n", "\t45 13519.866805 13495.866176 0.238079 0.240783\n", "\t50 13506.898673 13555.871313 0.219233 0.216347\n", "\t55 13549.458452 13571.589585 0.199181 0.197516\n", "\t60 13509.018362 13533.420573 0.185527 0.184918\n", "\t65 13486.072442 13516.614447 0.172556 0.171709\n", "\t70 13525.678046 13524.108128 0.157214 0.158169\n", "\t75 13482.036613 13504.251745 0.146788 0.146163\n", "\t80 13503.408657 13508.314541 0.135660 0.135793\n", "\t85 13518.750962 13538.331114 0.125845 0.125334\n", "\t90 13448.183311 13460.652047 0.118639 0.118766\n", "\t95 13466.878706 13463.743682 0.110304 0.111095\n", "\t100 13393.365097 13382.045279 0.105151 0.105867\n", "\t105 13320.206391 13378.572343 0.100640 0.099143\n", "\t110 13296.148156 13379.622903 0.095177 0.092730\n", "\t115 13277.749438 13274.010882 0.089847 0.089717\n", "\t120 13255.594266 13289.780395 0.084687 0.083850\n", "\t125 13216.769722 13209.527751 0.080475 0.080657\n", "\t130 13134.280941 13174.722631 0.077466 0.076697\n", "\t135 13143.270696 13135.588395 0.073214 0.073239\n", "\t140 13118.062526 13089.038929 0.069721 0.070019\n", "\t145 13129.778037 13107.025066 0.066201 0.066344\n", "\t150 13016.758140 13017.046389 0.064372 0.064109\n", "\t155 12952.608913 12909.922348 0.062374 0.062857\n", "\t160 12972.195659 12948.664891 0.058840 0.059189\n", "\t165 12901.782640 12903.756393 0.057205 0.057100\n", "\t170 12836.110004 12845.378254 0.055623 0.055307\n", "\t175 12813.578279 12803.645892 0.053283 0.053442\n", "\t180 12746.183324 12723.611287 0.052010 0.052088\n", "\t185 12642.400168 12634.318515 0.051181 0.051017\n", "\t190 12623.518838 12601.606254 0.049208 0.049256\n", "\t195 12545.983504 12565.203662 0.048173 0.047612\n", "\t200 12486.360144 12494.861637 0.046895 0.046639\n", "\t205 12395.242697 12397.928412 0.046184 0.045991\n", "\t210 12343.814221 12336.477565 0.045038 0.045007\n", "\t215 12244.749981 12249.970878 0.044401 0.044239\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 29/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21179e+43 erg / s Luminosity absorbed = 8.22822e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Plasma stratification:\n", "\t t_rad next_t_rad w next_w\n", "\tShell \n", "\t0 13233.150013 13208.310381 0.650204 0.657074\n", "\t5 13246.153944 13242.926346 0.561930 0.561964\n", "\t10 13321.909830 13330.763628 0.489863 0.489069\n", "\t15 13364.961838 13368.288133 0.435686 0.433913\n", "\t20 13391.543417 13384.044346 0.390502 0.389619\n", "\t25 13380.970341 13398.903276 0.354792 0.351207\n", "\t30 13424.798756 13448.522939 0.319017 0.316476\n", "\t35 13461.928210 13467.082069 0.291085 0.290285\n", "\t40 13503.499504 13508.888789 0.261981 0.263058\n", "\t45 13495.866176 13525.001658 0.240783 0.239508\n", "\t50 13555.871313 13562.425891 0.216347 0.217229\n", "\t55 13571.589585 13577.718177 0.197516 0.197810\n", "\t60 13533.420573 13597.504461 0.184918 0.181660\n", "\t65 13516.614447 13601.489270 0.171709 0.168048\n", "\t70 13524.108128 13568.698227 0.158169 0.156417\n", "\t75 13504.251745 13546.834728 0.146163 0.145109\n", "\t80 13508.314541 13555.479137 0.135793 0.133817\n", "\t85 13538.331114 13571.340270 0.125334 0.124472\n", "\t90 13460.652047 13472.928583 0.118766 0.119160\n", "\t95 13463.743682 13457.100412 0.111095 0.111430\n", "\t100 13382.045279 13428.556594 0.105867 0.105074\n", "\t105 13378.572343 13366.311341 0.099143 0.100007\n", "\t110 13379.622903 13379.869001 0.092730 0.093092\n", "\t115 13274.010882 13275.411938 0.089717 0.090079\n", "\t120 13289.780395 13246.539043 0.083850 0.085249\n", "\t125 13209.527751 13249.121435 0.080657 0.079777\n", "\t130 13174.722631 13135.663120 0.076697 0.077560\n", "\t135 13135.588395 13110.646608 0.073239 0.073567\n", "\t140 13089.038929 13087.479523 0.070019 0.070224\n", "\t145 13107.025066 13044.884829 0.066344 0.067281\n", "\t150 13017.046389 13037.484772 0.064109 0.063956\n", "\t155 12909.922348 12941.347490 0.062857 0.062526\n", "\t160 12948.664891 12942.316969 0.059189 0.059498\n", "\t165 12903.756393 12933.641487 0.057100 0.056895\n", "\t170 12845.378254 12859.107981 0.055307 0.055298\n", "\t175 12803.645892 12800.474636 0.053442 0.053562\n", "\t180 12723.611287 12750.762210 0.052088 0.052045\n", "\t185 12634.318515 12681.428533 0.051017 0.050617\n", "\t190 12601.606254 12600.917221 0.049256 0.049396\n", "\t195 12565.203662 12543.499820 0.047612 0.047945\n", "\t200 12494.861637 12490.808295 0.046639 0.046748\n", "\t205 12397.928412 12415.000297 0.045991 0.045719\n", "\t210 12336.477565 12336.243932 0.045007 0.044929\n", "\t215 12249.970878 12253.441375 0.044239 0.044120\n", "\n", " (\u001b[1mbase.py\u001b[0m:348)\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] t_inner 13000.000 K -- next t_inner 13000.000 K (\u001b[1mbase.py\u001b[0m:350)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Starting iteration 30/30 (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Luminosity emitted = 1.21579e+43 erg / s Luminosity absorbed = 8.19099e+42 erg / s Luminosity requested = 1.05928e+43 erg / s (\u001b[1mbase.py\u001b[0m:357)\n", "[\u001b[1mtardis.simulation.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Simulation finished in 30 iterations and took 302.03 s (\u001b[1mbase.py\u001b[0m:306)\n" ] } ], "source": [ "sim = run_tardis('blondin_model_compare_01.yml')" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/html": [ "Download Spectrum" ], "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "spectrum = pd.DataFrame(data = sim.transport.spectrum_integrated.frequency.value, columns=['frequency'])\n", "spectrum['luminosity'] = sim.transport.spectrum_integrated.luminosity_density_nu.value\n", "\n", "create_df_download_link(spectrum, title='Download Spectrum', filename='spectrum.csv', type='csv')" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "spec_artis_toy = pd.read_csv('spec_artis_toy.csv')" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/wkerzend/miniconda/envs/tardis3/lib/python3.6/site-packages/matplotlib/style/core.py:167: UserWarning: In /Users/wkerzend/.matplotlib/stylelib/talk_default.mplstyle: \n", "The text.latex.unicode rcparam was deprecated in Matplotlib 2.2 and will be removed in 3.1.\n", " styles = read_style_directory(stylelib_path)\n", " (\u001b[1mwarnings.py\u001b[0m:99)\n", "Populating the interactive namespace from numpy and matplotlib\n" ] }, { "data": { "application/javascript": "/* Put everything inside the global mpl namespace */\nwindow.mpl = {};\n\n\nmpl.get_websocket_type = function() {\n if (typeof(WebSocket) !== 'undefined') {\n return WebSocket;\n } else if (typeof(MozWebSocket) !== 'undefined') {\n return MozWebSocket;\n } else {\n alert('Your browser does not have WebSocket support.' +\n 'Please try Chrome, Safari or Firefox ≥ 6. ' +\n 'Firefox 4 and 5 are also supported but you ' +\n 'have to enable WebSockets in about:config.');\n };\n}\n\nmpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n this.id = figure_id;\n\n this.ws = websocket;\n\n this.supports_binary = (this.ws.binaryType != undefined);\n\n if (!this.supports_binary) {\n var warnings = document.getElementById(\"mpl-warnings\");\n if (warnings) {\n warnings.style.display = 'block';\n warnings.textContent = (\n \"This browser does not support binary websocket messages. \" +\n \"Performance may be slow.\");\n }\n }\n\n this.imageObj = new Image();\n\n this.context = undefined;\n this.message = undefined;\n this.canvas = undefined;\n this.rubberband_canvas = undefined;\n this.rubberband_context = undefined;\n this.format_dropdown = undefined;\n\n this.image_mode = 'full';\n\n this.root = $('
');\n this._root_extra_style(this.root)\n this.root.attr('style', 'display: inline-block');\n\n $(parent_element).append(this.root);\n\n this._init_header(this);\n this._init_canvas(this);\n this._init_toolbar(this);\n\n var fig = this;\n\n this.waiting = false;\n\n this.ws.onopen = function () {\n fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n fig.send_message(\"send_image_mode\", {});\n if (mpl.ratio != 1) {\n fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n }\n fig.send_message(\"refresh\", {});\n }\n\n this.imageObj.onload = function() {\n if (fig.image_mode == 'full') {\n // Full images could contain transparency (where diff images\n // almost always do), so we need to clear the canvas so that\n // there is no ghosting.\n fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n }\n fig.context.drawImage(fig.imageObj, 0, 0);\n };\n\n this.imageObj.onunload = function() {\n fig.ws.close();\n }\n\n this.ws.onmessage = this._make_on_message_function(this);\n\n this.ondownload = ondownload;\n}\n\nmpl.figure.prototype._init_header = function() {\n var titlebar = $(\n '
');\n var titletext = $(\n '
');\n titlebar.append(titletext)\n this.root.append(titlebar);\n this.header = titletext[0];\n}\n\n\n\nmpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n\n}\n\n\nmpl.figure.prototype._root_extra_style = function(canvas_div) {\n\n}\n\nmpl.figure.prototype._init_canvas = function() {\n var fig = this;\n\n var canvas_div = $('
');\n\n canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n\n function canvas_keyboard_event(event) {\n return fig.key_event(event, event['data']);\n }\n\n canvas_div.keydown('key_press', canvas_keyboard_event);\n canvas_div.keyup('key_release', canvas_keyboard_event);\n this.canvas_div = canvas_div\n this._canvas_extra_style(canvas_div)\n this.root.append(canvas_div);\n\n var canvas = $('');\n canvas.addClass('mpl-canvas');\n canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n\n this.canvas = canvas[0];\n this.context = canvas[0].getContext(\"2d\");\n\n var backingStore = this.context.backingStorePixelRatio ||\n\tthis.context.webkitBackingStorePixelRatio ||\n\tthis.context.mozBackingStorePixelRatio ||\n\tthis.context.msBackingStorePixelRatio ||\n\tthis.context.oBackingStorePixelRatio ||\n\tthis.context.backingStorePixelRatio || 1;\n\n mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n\n var rubberband = $('');\n rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n\n var pass_mouse_events = true;\n\n canvas_div.resizable({\n start: function(event, ui) {\n pass_mouse_events = false;\n },\n resize: function(event, ui) {\n fig.request_resize(ui.size.width, ui.size.height);\n },\n stop: function(event, ui) {\n pass_mouse_events = true;\n fig.request_resize(ui.size.width, ui.size.height);\n },\n });\n\n function mouse_event_fn(event) {\n if (pass_mouse_events)\n return fig.mouse_event(event, event['data']);\n }\n\n rubberband.mousedown('button_press', mouse_event_fn);\n rubberband.mouseup('button_release', mouse_event_fn);\n // Throttle sequential mouse events to 1 every 20ms.\n rubberband.mousemove('motion_notify', mouse_event_fn);\n\n rubberband.mouseenter('figure_enter', mouse_event_fn);\n rubberband.mouseleave('figure_leave', mouse_event_fn);\n\n canvas_div.on(\"wheel\", function (event) {\n event = event.originalEvent;\n event['data'] = 'scroll'\n if (event.deltaY < 0) {\n event.step = 1;\n } else {\n event.step = -1;\n }\n mouse_event_fn(event);\n });\n\n canvas_div.append(canvas);\n canvas_div.append(rubberband);\n\n this.rubberband = rubberband;\n this.rubberband_canvas = rubberband[0];\n this.rubberband_context = rubberband[0].getContext(\"2d\");\n this.rubberband_context.strokeStyle = \"#000000\";\n\n this._resize_canvas = function(width, height) {\n // Keep the size of the canvas, canvas container, and rubber band\n // canvas in synch.\n canvas_div.css('width', width)\n canvas_div.css('height', height)\n\n canvas.attr('width', width * mpl.ratio);\n canvas.attr('height', height * mpl.ratio);\n canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n\n rubberband.attr('width', width);\n rubberband.attr('height', height);\n }\n\n // Set the figure to an initial 600x600px, this will subsequently be updated\n // upon first draw.\n this._resize_canvas(600, 600);\n\n // Disable right mouse context menu.\n $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n return false;\n });\n\n function set_focus () {\n canvas.focus();\n canvas_div.focus();\n }\n\n window.setTimeout(set_focus, 100);\n}\n\nmpl.figure.prototype._init_toolbar = function() {\n var fig = this;\n\n var nav_element = $('
')\n nav_element.attr('style', 'width: 100%');\n this.root.append(nav_element);\n\n // Define a callback function for later on.\n function toolbar_event(event) {\n return fig.toolbar_button_onclick(event['data']);\n }\n function toolbar_mouse_event(event) {\n return fig.toolbar_button_onmouseover(event['data']);\n }\n\n for(var toolbar_ind in mpl.toolbar_items) {\n var name = mpl.toolbar_items[toolbar_ind][0];\n var tooltip = mpl.toolbar_items[toolbar_ind][1];\n var image = mpl.toolbar_items[toolbar_ind][2];\n var method_name = mpl.toolbar_items[toolbar_ind][3];\n\n if (!name) {\n // put a spacer in here.\n continue;\n }\n var button = $('